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1. 


PAQCAGE  ORGANIZATION 


1.1  FLOW  CHART 

The  operation  of  REXOR  in  the  broadest  sense  consists  of  first  establishing 
an  equilibrium  set  of  initial  conditions  followed  by  a time  history  se- 
quence determined  by  control  inputs  and  the  equations  of  mcticn.  The  first 
part  of  the  operation  is  called  TRIM  and  the  second  part  is  termed  FLY. 

The  program  operations  REXOR  computes  in  TRIM  are  given  in  Figure  1-1.  The 
FLY  operations  are  shown  in  Figure  1-2.  The  computa'ion  blocks  are  briefly 
annotated  for  completeness  in  these  figures,  and  are  developed  more  fully 
in  the  following  subsections.  Likewise  the  available  options  are  explained 
below. 


1.2  C.PERATION  MODES 
1.2.1  TRIM 

Referring  to  Figure  1-1,  the  first  major  operation  performed  in  REXOR  is 
for  the  executive  routine  MAIN  to  call  the  subroutine  TRIM.  As  prograimed, 
TRIM  can  process  the  equations  of  motion  of  the  main  rotor  ensemble  and 
fuselage  accelerations.  The  main  rotor  trim  configuration  is  the  same  as 
will  be  used  in  FLY,  that  is,  a complete  rctor  and  control  system.  The 
code  in  TRIM  contains  the  equations  for  other  options.  However,  these  are 
not  fully  developed  or  sufficiently  reliable  to  be  considered  as 
operational. 

The  subroutines  LOADS,  SWEEP,  ACCEL  and  associated  subroutines  are  called 
from  TRIM.  These  subroutines  form  the  generalized  mass  and  force  matrices 
as  developed  in  Volume  I,  Section  6 for  the  rotor  and  control  system.  The 
correcting  acceleration.  Section  6.3,  is  formed  by  the  subroutine  MIC26. 

The  predictor  numericeil  integration  is  done  in  the  subroutine  INTG  with  the 
entry  points  of  PRED  and  CORECT. 

The  subroutine  TRIM  uses  the  newly  updated  set  of  accelerations  to  adjust 
the  selected  control  .set  to  converge  on  the  desired  initial  condition  set. 
The  initial  condition  set  can  specify  uns,ccelerated  flight  or  an  initial 
load  factor  in  a coordinated  turn.  "Hie  control  set  selection  available  is 
given  in  Volume  III.  The  trim  error  acceleration  set  (determined  by  the 
control  set  selected)  operates  the  control  set  through  input  gain  factors 
to  null  the  errors.  Trim  is  detected  by  simultaneo  steadiness  of  all  the 
selected  trim  controls.  The  exit  from  TRIM  is  to  MAIN  which  in  turn  trans- 
fers control  directly  to  FLY.  A control  flag  to  save  the  completed  trim 
data  is  available  and  operates  prior  to  the  transfer  to  MAIN  and  FLY.  Tlie 
feature  saves  the  converge.uce  time  when  running  the  same  case  at  a later 
date. 
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1.2.2  FLY 


As  shown  in  Figure  1-2,  the  operation  of  FLY  closely  resembles  that  of 
TRIM;  namely  the  equations  of  motion  are  still  formed  by  the  subroutines 
LOADS,  SWEEP,  ACCEL  and  associated  routines.  The  difference  is  that  the 
entire  set  of  equations  ore  being  coraputed  and  that  the  controls  are 
driven  by  commands  rather  than  trim  error  balemce  sources. 

The  control  system  can  operate  with  a number  of  different  configurations. 

A hard  swashplate  and  flexible  swashplate  - external  control  gyro  configu- 
rations are  computed  directly  in  the  subroutine  FLY.  For  the  isolated  in- 
ternal gyro  system  (Lockheed  Advanced  Mechanical  Control  System)  the  sub- 
routine IGYRO  is  used.  This  subroutine  in  turn  uses  the  subroutine  BIRD! 
(with  multiple  entry  points)  to  form  and  solve  the  equations  of  motion  of 
this  system  apart  from  the  main  flow  of  computation.  The  shaft  bending 
equations  are  aisj  integrated  by  this  branch  computation. 

The  input  commands  for  cyclic,  collective  and  pedals  come  from  the  sub- 
routine CNTROL  (entry  point  PICIfTL).  This  gives  a tine  history  input  posi- 
tion and  rate  commands  from  the  input  data  set  (Volume  III).  A number  of 
flight  profile  follower  - autopilots  exist  in  the  code,  and  can  be  used  at 
this  point.  However,  the  code  is  not  considered  to  be  documented  or 
operational. 

The  available  airfoil  data  schemes  are  shown  in  Figure  1-2.  For  clarity 
the  less  involved  7 table  lookup  is  shown  in  the  computation  loop.  The 
fast  aero  lookup  system  is  shown  as  an  alternate  at  the  bottom  of  the  fig- 
ure together  with  the  dynamic  stall  option  (subroutine  STALL)  which  may 
be  used  with  this  lookup  procedure.  The  background  material  for  the  lookup 
routines  is  contained  in  Volume  I,  Section  7 and  Section  b of  this  Volume. 

Other  options  which  are  activated  from  the  subroutine  fly  are  variable 
rotor  speed  (subroutine  ETORQ),  quasi-static  blade  torsion  (subroutine 
TORSI  entry  TORS),  and  isolated  pitch  horn  bending  (subroutine  PHORN). 

Operation  of  the  subrout*  ne  FLY  continues  in  a loop  fashion  until  the 
allowed  maneuver  time  limit  (input)  or  unrealistic  loads  are  detected  iue 
to  true  rotorcraft  problems  or  computation  numerical  difficulties.  Exit 
from  FLY  returns  control  to  MAIN.  Here,  either  a new  case  or  an  execution 
termination  occurs  depending  on  the  input  data. 


1. 3 ROUTINES  IN  PACKAGE 

In  this  section  the  computer  routines  which  comprise  REXOR  are  explained  in 
sufficient  detail  to  allow  the  reader  to  identify  and  map  out  the  computa- 
tion procedures  explained  in  Volume  I.  The  REXOR  routines  sure  summarized 
in  Table  1-1.  The  following  subheadings  treat  the  routines  as  listed  in 
this  table.  The  subroutines  ACCEL,  LOADS  and  SWEEP  form  the  nucleus  of  the 
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TABLE  1-1.  REXOR  ROUTINES 

Routine 

Name 

Entry  Points 

1- 

ACCELl 

ACCELl,  ACCEL 

2, 

AERO 

AERO 

3. 

BIRDl 

BIRDl,  PBIRD,  ABIRD,  BBIRD,  COR 

1*. 

BMOVE 

BMOVE, 

5. 

BSCALE 

BSCALE 

6. 

CMLOOK 

CMLOOK 

7. 

CNTROL 

CNTROL,  APCKTL,  PICNTL 

8. 

COULOM 

COULOM 

9. 

DERTV 

DERIV,  DERIVE 

10. 

DWASH 

DV/ASH 

11. 

ETORQl 

ETORQl,  ETORQ 

12. 

FLY 

FLY 

13. 

HAR141 

HARl-0.,  HARM,  HARM2 

li*. 

IGYRO 

IGYRO 

15. 

INTO 

INTG,  PRED,  COKECT 

l6. 

LOADSl 

LOADSl,  LOADS 

17. 

LSTDAT 

LSTDAT 

13. 

MAIN 

(CONTROLLING  ROUTINE) 

19. 

MIC26 

MIC26 

20. 

MINVR 

MINVR 

21. 

MPRNT 

f-lPRNT 

22. 

PDATE 

PDATE 

23. 

PHORN 

PHORN 

2k. 

PILOTA 

PILOTA,  PTLOTI 

25. 

PRINTl 

PRINTl,  PRINT 

26. 

PROP 

(BLOCK  DATA) 

27. 

RCPLOT 

RCPLOT 

28. 

RE  ADI  N 

READIN 

29. 

STALL 

STALL 
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Preceiiing  page  blank 


TABLE  1-1  - Continued 


Routine 


Name 

Entry  Points 

30. 

SWEEPl 

SWEEPl,  SWEEP 

31. 

TORSI 

TORSI,  TORS 

32. 

TRIM 

TRIM 

33. 

TRMPUn 

TRMPUN,  DSHIFT 

3h. 

XTERMO 

XTERMO,  XTERMl,  XTERM2,  XTERM3,  XTERM1+,  XTERM5, 

XTERM6,  XTERM7,  XTERM8,  XTERM9,  XTERMA,  XT'ERKC, 
XTERMD,  XTERME 


3i). 

XTRPl 

XTRPl 

36. 

XTRP2 

XTRP2 

37. 

XTRP3 

XTRP3 

38. 

CMA 

(BLOCK  DATA) 

39. 

CMAM 

(BLOCK  DATA) 

ho. 

CLAT 

(BLOCK  DATA) 

1»1. 

CLATM 

(BLOCK  DATA) 

h?.. 

CDAT 

(BLOCK  DATA) 

J*3. 

CDATM 

(BLOCK  DATA) 

hh. 

XTRP 

XTRP 

h^. 

XTRPl* 

XTRPl* 

hS, 

XTRPCA 

XTRPCA 

hi. 

CMTICL 

(BLOCK  DATA) 

SIMILAR  BLOCKS  -2  TO 

1*8. 

CMTICD 

(BLOCK  DATA) 

h9. 

CTISCL 

(BLOCK  DATA) 

SIMILAR  BLOCKS  -2  TO  -7 

50. 

CTISCD 

(BLOCK  DATA) 

51. 

CMTICM 

(BLOCK  DATA) 

52. 

CMT2CM 

(BLOCK  DATA) 
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computation  procedure  known  as  REXOR.  Therefore  there  subroutines  are 
covered  in  detail,  and  matched  to  the  symbology  and  equations  of  Volume  I. 
Another  class  of  subroutines  contain  computation  subgroupings  required  by 
the  routines  just  mentioned,  or  setup  the  integration  and  equation  solving 
process.  These  routines  are  coveied  in  less  detail,  but  with  enough  ex- 
planatory material  for  the  reader  to  see  how  the  programming  accomplishes 
the  desired  task.  The  last  class  of  routines  is  one  of  utility  functions 
s ich  as  plotting,  printing  and  date.  These  entries  are  noted  as  to  input- 
0".tput  function. 

1.3.1  ACCELl 


As  shown  in  Table  1-1,  the  subroutine  ACCELl  follows  the  pattern  of  many 
REXOR  routines  in  that  the  first  entry  point  initializes  often  used  con- 
stants, and  the  subsequent  entries  perfom  the  required  computations.  The 
entry  point  ACCELl  sets  the  operating  RA  equivalence.  Numerous  sines  and 
cosines  are  precomputed  and  saved  as  well  as  conversion  of  degree  data  to 
radians.  Many  diametral  inertias  are  computed  as  half  the  data  entry  polar 
inertia  values. 

The  entry  point  ACCEL  in  conjunction  with  similar  points  LOADS  and  SWEEP 
form  the  computation  nucleus  of  REXOR.  ACCEL  gathers  the  information  to 
form  the  generalized  mass  and  force  matricies,  and  controls  the  accelera- 
tion update  sequence.  The  majority  of  the  development  of  Volume  I,  Sec- 
tion 6,  except  for  blade  geometry  and  summation,  is  coded  in  this  entry. 

As  shown  in  Figure  1-3,  the  entry  ACCEL  collects  the  generalized  masses  and 
forces  to  proceed  with  determining  acceleration  terms.  When  REXOR  is  in 
the  TRIM  mode  of  operation,  the  mass  and  force  equations  are  solved  twice 
per  pass  in  the  subroutine.  First,  only  the  blade  degree  of  freedom  equa- 
tions are  solved  for  accelerations  to  be  used  in  the  time  integration. 
Second,  the  entire  degree  of  freedom  acceleration  vector  is  solved.  The 
TRIM  acceleration  terms  are  used  from  this  second  trial  to  readjust  the 
controls  to  iterate  a stable  trimmed  flight  condition.  Two  tria.  3 are 
needed  so  that  the  blade  accelerations  used  for  blade  integrations  are  not 
contaminated  with  incorrect  accelerations  from  the  other  degrees  of  freedom. 
In  the  FLY  mode  the  entire  acceleration  vector  is  developed  and  integrated. 
Therefore,  the  solution  sequence  takes  place  only  once  per  pass  in  ACCEL. 

The  generalized  masses  ai*e  coded  as  QMG  and  the  generalized  forces  eire 
coded  as  the  QFG  array.  The  notation  of  these  arrays  is  locked  to  the 
problem  variable  names  as  given  in  Table  1-2.  The  coding  may  be  traced 
back  to  Volume  I by  use  of  this  table,  keeping  in  mind  that  the  mass  matrix 
is  symmetrical.  For  example  QMG(21,  l6)  is  M^  „ from  Volume  I, 

Section  6.9.  h°PH'* 

The  generalized  masses  are  assembled  from  precomputed  inertias,  partial 
derivatives,  subcomputations  and  results  from  SWEEP.  The  inertias,  center 
of  gravity  offsets,  and  masses  can  be  directly  read  from  the  code  noting 


Figure  1-3-  ACCEL  Computation  Flow. 
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TABLE  1-2.  QMG  AM)  QFG  AHRAY  CODING 


Degree  of  Freedom 


Blade  1 


Blade  2 


Blade  3 


Blade  k 


Swash  Plate 


Rotor  Azimuth 


Principal  Ref.  Axis 


Shaft  Bending 


Gyro  Variables  4>q,  6q  Are  Solved  For  Independently 
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that  F stands  for  fuselage  and  B stands  for  princi"-ul  reference  axis 
(H  for  hub).  The  B notation  stems  from  the  old  refe.rr-nce  B for  Body.  BAR 
means  (“),  and  identifies  a center-of-gravity  offset. 

Partial  derivatives  are  grouped,  and  use  a six-letter  identifier  generally 
starting  with  P for  pai’tial.  The  sequence  of  letters  gives  the  partial  de- 
rivative numerator  and  denominator.  For  example  PZFPYB  is  partial  Z fuse- 
lage, with  respect  to  partial  Y body  or  ^Zp/3Yf^.  The  six-letter  limit 
calls  for  improvisation.  PHFPSB  is  and  THFTHB  is  36p/30{|.  Terms 

involving  blade  partial  derivatives  which  require  spanwise  blade  integra- 
tion are  transferred  from  SWEEP  in  the  array  F which  is  described  under  the 
heading  SWEEPl.  A derivative  listing  is  given  in  Table  1-3. 

A number  of  subcomputations  grouped  as  BM and  CM are  made  for  body 

(hub)  mass  and  control  mass  respectively.  The  mnemonics  used  are  not 
directly  readable,  but  should  be  referenced  to  the  QMG  calculations  for 
meaning. 

The  QFG  array  is  assembled  from  blade  data  (also  in  the  F array  mentioned), 
data  from  LOADS  and  swashplate  loads.  The  latter  are  developed  within 
ACCEL  fr<m  Volume  I,  Section  6.10.  The  array  QLOADS  from  LOADS  is  the  non- 
main  rotor  aerodynamic  loads  as  developed  in  Volume  I,  Section  7.U.  The 
QFG  components  involving  these  loads  are  assembled  with  the  appropriate 
inertias,  partials,  etc.,  as  the  equations  require. 

Using  the  con5)leted  QMG  and  QFG  arrays,  a new  correction  acceleration 
vector  DELA  is  found  from  the  subroutine  MIC26.  DELA  then  is  added  to  the 
running  estimate  acceleration  vector  YDD.  The  YDD  array  has  the  same 
ordering  as  given  in  Table  1-2. 

Due  to  the  method  of  introducing  the  gravity  vector,  discussed  in  Volume  I, 
Section  5. 5.1,  gravity  acceleration  as  well  as  maneuvering  loads  will 
appear  on  the  vector  triad  YDD(21,  22,  23).  These  accelerations  are 
balanced  by  the  force  vector,  thus  propagating  the  acceleration  due  to 
gravity  throughout  the  problem.  These  accelerations  do  not  need  to  be  in- 
tegrated for  use  elsewhere,  therefore  the  positions  YD(21,  22,  23)  are 
used  for  the  corresponding  accelerations  without  gravity.  The  multiple  use 
of  the  YDD,  YD,  Y array  is  not  limited  to  the  gravity  terras,  an  array  sum- 
mary is  given  in  Table  l-i<  for  dimensions  greater  than  20,  and  will  be 
referred  to  under  subsequent  headings.  No  foldover  exists  for  dimensions 
less  than  21,  and  Table  1-2  may  be  referenced, 

1.3.2  AERO 

The  subroutine  AERO  is  called  in  the  Fast  Aero  main  rotor  aerodynamic  data 
lookup  procedure.  The  operations  of  this  procedure  are  developed  in  Vol- 
ume I,  Section  7.2.U’,  Section  k of  this  volume; and  Volume  III,  Sec- 
tion 3. 3.7. 3.  This  Fast  Aero  procedure  uses  airfoil  data  which  has  been 
interpolated  into  equally  spaced  increments,  and  reduced  only  to  the  types 
and  sections  of  blade  data  actually  needed.  The  set  currently  in  REXOR  is 
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TABLE  1-3.  PARTIAL  DERIVATIVES  IN  ACCEL 


Shaft  Bending 

Symbology  FORTRAN  Name 


r 

* /\n 

- 

OF 

3d) 

► = ^ 

PYFPPS  = YPHIF  = RA(365) 

S 

- 

J 

^ ■« 
XTHTF  = RA(36U) 


V 


’ - 

PHFPKS 
< _ • 

PSIPKS 

>■ 

PHFTHS 
< THFTHS  ' 
PSITHS 


Svashplate  Motions 


(1-2) 


(1-3) 


(1-1*) 


• 

PrCiOF,  PFGGD,  - 

^^SP 

PTGGK,  PTGGD, 

9,. 

L ^sp . 

. PSIGGK,  PSIGGD.  - . 

3z 


OSP 


^SP 


1 


(1-5) 

(1-6) 


TABLE  1-3  - Continued 


Swashplate  Motions  (Continued) 

Symbology  FORTRAN  Neime 


3 (p 


FN 


SP 


3 p 


FN 


'SP 


PFGK(I) 


PFGD(I) 


(1-7) 


(1-8) 


3 ({i 


FN 


7 

-SP 


PFHG(I) 


Pitch  Horn  Bending 


3({) 

^FN 


PPEPH 


PHn 


Principal  Axis  Motions 


(1-9) 


(1-10) 


3v. 

’ PXRXB(I),  PXRYBd),  - 

OBLn 
S r 

= 

PYRXB(I),  PYRYBd), 

- , - , PZRZB(I). 

(1-11) 


3,. 

■PXFPXS,  PXFPYB,  PXFPZB' 

OF 

= 

PYFPXB,  PYFPYB,  PYFPZB 

^ r 

U H J 

, PZFPXB,  PZFPYB,  PZFPZB. 

(1-12) 
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TABLE  1-1.  REXOR  ROUTINES 

Routine 

Name 

Entry  Points 

1. 

ACCELl 

ACCELl,  ACCEL 

2. 

AERO 

AERO 

3. 

BIRDl 

BIRDl,  PBIRD,  ABIRD,  BBIRD,  COR 

h. 

BMOVE 

BMOVE, 

5. 

BSCALE 

BSCALE 

6. 

CMLOOK 

CMLOOK 

7. 

CNTROL 

CNTROL,  APCNTL,  PICNTL 

8. 

COULOM 

COULOM 

9. 

DERIV 

DERIV,  DERIVF 

10. 

DWASH 

DWASH 

11. 

ETORQl 

ETORQl,  ETORQ 

12. 

FLY 

FLY 

13. 

HARMl 

HARMl,  HARM,  HARM2 

lU. 

IGYRO 

IGYRO 

15. 

INTO 

INTO,  PRED,  CORECT 

16. 

LOADSl 

LOADSl,  LOADS 

17. 

LSTDAT 

LSTDAT 

18. 

MAIN 

(CONTROLLING  ROUTINE) 

19. 

1CEC26 

MIC26 

20. 

MINVR 

MINVR 

21. 

MPRNT 

MPRNT 

22. 

PDATE 

PDATE 

23. 

PHORN 

PHORN 

2k. 

PILOTA 

PILOTA,  PILOT! 

25. 

PRINTl 

PRINTl,  PRINT 

26, 

PROP 

(BLOCK  DATA) 

27. 

RCPLOT 

RCPLOT 

28. 

READIN 

READIN 

29. 

STALL 

STALL 
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TABLE  1-1  - Continued 


Routine 


Name 

Entry  Points 

30.  SWEEPl 

SWEEPl,  S\7EEP 

31.  TORSI 

TORSI,  TORS 

32.  TRIM 

TRIM 

33.  TRMPUH 

TRMPUN,  DSHIFT 

3i*.  XTERMO 

:(TERM0,  XTERMl,  XTERM2,  XTERM3,  XTERM4,  XTERM5 
XTERM6,  XTERM7,  XTERM8,  XTERM9,  XTERMA,  XTERMC 

XTERMD,  XTERME 

35.  XTRPl 

XTRPl 

36.  XTRP2 

XTRP2 

37.  XTRP3 

XTRP  3 

38.  CMA 

(BLOCK  DATA) 

39.  CMAM 

(BLOCK  DATA) 

i*0.  CLAT 

(BLOCK  DATA) 

1*1.  CLATM 

(BLOCK  DATA) 

1*2.  CDAT 

(BLOCK  DATA) 

1*3.  CDATM 

(BLOCK  DATA) 

1*1*.  XTRP 

XTRP 

1*5.  XTRPl* 

XTRP4 

46.  XTRPCA 

XTRPCA 

47.  CMTICL 

(BLOCK  DATA) 

SIMILAR  BLOCKS  -2  TO  -7 

48.  CMTICD 

(BLOCK  DATA) 

49.  Cl'ISCL 

(BLOCK  DATA) 

SIMILAR  BLOCKS  -2  TO  -7 

50.  CTlu^c 

(BLOCK  DATA) 

51.  CMTICM 

(BLOCK  DATA) 

52.  CMT2CM 

(BLOCK  DATA) 

r 

K 

»& 

k 


computation  procedure  known  as  REXOR.  Therefore  these  subroutines  are 
covered  in  detail,  and  matched  to  the  symbology  and  equations  of  Volume  I. 
Another  class  of  subroutines  contain  computation  subgroupings  required  by 
the  routines  Just  mentioned,  or  setup  the  integration  and  equation  solving 
process.  These  routines  are  covered  in  less  detail,  but  with  enough  ex- 
planatory material  for  reader  to  see  how  the  programming  accomplishes 
the  desired  task.  The  last  class  of  routines  is  one  of  utility  functions 
such  as  plotting,  printing  and  date.  These  entries  are  noted  as  to  input- 
output  function. 

1.3.1  ACCELl 


As  shown  in  Table  1-1,  the  subroutine  ACCELl  follows  the  pattern  of  many 
REXOR  routines  in  that  the  first  entry  point  initializes  often  used  con- 
stants, and  the  subsequent  entries  perform  the  required  computations.  The 
entry  point  ACCELl  sets  the  operating  RA  equivalence.  Numerous  sines  and 
cosines  are  precomputed  and  saved  as  well  as  conversion  of  degree  data  to 
radians.  Many  diametral  i.iertias  are  computed  as  half  the  data  entry  polar 
inertia  values. 

The  entry  point  ACCEL  in  conjunction  with  similar  points  LOADS  and  SWEEP 
form  the  computation  nucleus  of  REXOR.  ACCEL  gathers  the  information  to 
form  the  generalized  mass  and  force  matricies,  and  controls  the  accelera- 
tion update  sequence.  The  majority  of  the  development  of  Volume  I,  Sec- 
tion 6,  except  for  blade  geometry  and  summation,  is  coded  in  this  entry. 

As  shown  in  Figure  1-3,  the  entry  ACCEL  collects  the  generalized  masses  and 
forces  to  proceed  with  determining  acceleration  terms.  When  REXOR  is  in 
the  TRIM  mode  of  operation,  the  mass  and  force  equations  are  solved  twice 
per  pass  in  the  subroutine.  First,  only  the  blade  degree  of  freedom  equa- 
tions are  solved  for  accelerations  to  be  used  in  the  time  integration. 
Second,  the  entire  degree  of  freedom  acceleration  vector  is  solved.  The 
TRIM  acceleration  terms  are  used  from  this  second  trial  to  readjust  the 
controls  to  iterate  a stable  trimmed  flight  condition.  Two  trials  are 
needed  so  that  the  blade  accelerations  used  for  blade  integrations  are  not 
contaminated  with  incorrect  accelerations  from  the  other  degrees  of  freedom. 
In  the  FLY  mode  the  entire  acceleration  vector  is  developed  and  integrated. 
Therefore,  the  solution  sequence  teikes  place  only  once  per  pass  in  ACCEL. 

The  generalized  masses  are  coded  as  QMG  and  the  generalized  forces  are 
coded  as  the  QFG  array.  The  notation  of  these  arrays  is  locked  to  the 
problem  variable  names  as  given  in  Table  1-2.  The  coding  may  be  traced 
back  to  Volume  I by  use  of  this  table,  keeping  in  mind  that  the  mass  matrix 
is  symmetrical.  For  example  QMG(21,  l6)  is  M^  „ from  Volume  I, 

Section  6.9.  PHt 

The  generalized  marses  are  assembled  from  precomputed  inertias,  partial 
derivatives,  subcoraputations  and  results  from  SWEEP.  The  inertias,  center 
of  gravity  offsets,  and  masses  can  be  directly  read  from  the  code  noting 
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Figure  1-3.  ACCEL  Computation  Flow. 
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TABLE  1-2.  QMS  AOT)  QFG  ASflAY  CODING 


Degree  of  Freedom 


Blade  1 


Blade  2 


Blade  3 


Blade  i* 


Swash  Plate 


Rotor  Azimuth 


Principal  Ref.  Axis 


Shaft  Bending 


Gyro  Variables 


0„  Are  Solved  For  Independently 
G 
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that  F stands  for  fuselage  and  B stands  for  principal  reference  axis 
(H  for  hub).  The  B notation  stems  from  the  old  reference  B for  Body.  BAB 
means  (~),  and  identifies  a center-of-gravity  offset. 

Partial  derivatives  are  grouped,  and  use  a six-letter  identifier  generally 
starting  with  P for  partial.  The  sequence  of  letters  gives  the  partial  de- 
rivative numerator  and  denominator.  For  example  PZFPYB  is  partial  Z fuse- 
lage, with  respect  to  partial  Y body  or  9Zp/3Yjj.  The  six-letter  limit 
calls  for  improvisation.  PHFPSB  is  3(j)p/3'}'jj  and  THFTHB  is  96p/99j£.  Terms 
involving  blade  partial  derivatives  which  require  spanwise  blade  integra- 
tion are  transferred  from  SWEEP  in  the  array  F which  is  described  under  the 
heading  SWEEP3 . A derivative  listing  is  given  in  Table  1-3. 

A number  of  subcomputations  grouped  as  BM and  CM are  made  for  body 

(hub)  mass  and  control  mass  respectively.  The  mnemonics  used  are  not 
directly  readable,  but  should  be  referenced  to  the  QMG  calculations  for 
meaning. 

The  QFG  array  is  assembled  from  blade  data  (also  in  the  F array  mentioned), 
data  from  LOADS  and  swashplate  loads.  The  latter  are  developed  within 
ACCEL  from  Volume  I,  Section  6.10.  The  array  QLOADS  from  LOADS  is  the  non- 
main  rotor  aerodynamic  loads  as  developed  in  Volume  I,  Section  7.^.  The 
QFG  components  involving  these  loads  are  assembled  with  the  appropriate 
inertias,  partials,  etc.,  as  the  equations  require. 

Using  the  completed  QMG  and  QFG  arrays,  a new  correction  acceleration 
vector  DELA  is  found  from  the  subroutine  MIC26.  DELA  then  is  added  to  the 
running  estimate  acceleration  vector  YDD.  The  YDD  array  has  the  same 
ordering  as  given  in  Table  1-2. 

Due  to  the  method  of  introducing  the  gravity  vector,  discussed  in  Volume  I, 
Section  5«5«1,  gravity  acceleration  as  well  as  maneuvering  loads  will 
appear  on  the  vector  triad  YDD(21,  22,  23).  These  accelerations  are 
balanced  by  the  force  vector,  thus  propagating  the  acceleration  due  to 
gravity  throughout  the  problem.  These  accelerations  do  not  need  to  be  in- 
tegrated for  use  elsewhere,  therefore  the  positions  YD(21,  22,  23)  are 
used  for  the  corresponding  accelerations  without  gravity.  The  multiple  use 
of  the  YDD,  YD,  Y array  is  not  limited  to  the  gravity  terms.  An  array  sum- 
mary is  given  in  Table  1-h  for  dimensions  greater  than  20,  and  will  be 
referred  to  under  subsequent  headings.  No  foldover  exists  for  dimensions 
less  than  21,  and  Table  1-2  may  be  referenced. 

1.3.2  AERO 

The  subroutine  AERO  is  called  in  the  Fast  Aero  main  rotor  aerodynaiiiic  data 
lookup  procedure.  The  operations  of  this  procedure  are  developed  in  Vol- 
ume I,  Section  7.2.h;  Section  !i  of  this  volume;  and  Volume  III,  Sec- 
cion  3.3.T.3.  This  Fast  Aero  procedure  uses  airfoil  data  which  has  been 
interpolated  into  equally  spaced  increments,  and  reduced  only  to  the  tj^jes 
and  sections  of  blade  data  actually  needed.  The  set  currently  in  REXOR  is 
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TABLE  1-3  - Continued 


Swashplate  Motions  (Continued) 

Symbology  FORTRAW  Name 


PFGKd^ 


(2-7) 


3 


0 


SP 


PFGD(I) 


(1-8) 


3 


^SP 


PFHG(I) 


(1-9) 


Pitch  Horn  Bending 


9b 

PKn 


PFBPH 


(1-10) 


Principal  Axis  Motions 


3r 

* PXRXB(I),  PXRYB(I), 

OBLn 
9 r 

= 

PYRXB(I),  PYRYB(I), 

- , - , PZRZB(I). 

' PXFPXB, 

PXFPYB, 

pxfpzb" 

OF 

PYFPXB, 

PYFPYB, 

PYFPZB 

^ r 

L H ^ 

, PZFPXB, 

PZFPYB, 

PZFPZB, 
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TABLE  1-3  - Continued 


Principal  Axis  Motions  (Continued) 


Symbology 


FORTRAN  Name 


- , PXFPTB,  PXFPSB 

PYFPPB,  - , PYFPSB 


PHFPHS»  PHFTHB,  PHFPSB 
- , THFTHB,  THFPSB 

PSFPHB,  PSFTHB,  PSFPSB 


PFGFB,  PFGTB,  PFGSIB 
- , PTGTB,  PTGSIB 

PSIGFB,  PSIGTB,  PSIGSB 


Rotor  (Engine)  Rotation 


PFSI(I) 


(1-16) 
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TABLE  1-U.  MOTION  ARRAY  FOLDOVER 


Array  Contents 

I'DD 

YD 

Y 

• • 

2* 

• 

w„ 

w 

H 

« 

H 

H 

% 

Ph 

Ph 

% 

% 

'Ih 

^H 

^H 

! = ZDD(3) 

<^E 

<j>E 

Og  = ZDD(U) 

®E 

- 

• 

• 

for  the  AH-56A  Cheyenne.  However  any  set  of  tables  can  be  inserted  so 
long  as  the  call  argument  and  table  arrangement  are  correct. 

AERO  references  the  blade  data  stored  in  block  data  sets  CLATM,  CDATM» 
CLAT,  CDAT,  CMAM,  CMA.  These  block  data  sets  are  entered  with  the  inter- 
polation routines  XTRPl,  XTRP2  and  XTRP3  depending  on  whether  1,  2 or  3 
parameters  are  being  interpolated. 

The  subroutine  call  argument  is  as  follows: 

Input: 

XN  Normalized  blade  span  location  from  0 to  1. 

MACH  Local  section  Mach  number. 

ALFAR  Section  angle  of  attack,  rad. 

ICMBL  Blade  pitching  moment  coefficient  flag.  This  element  is 
also  referred  to  as  IBLADE  in  the  calling  routines,  and 
sets  the  trailing-edge  configuration  for  AH-56A  blades. 

DCf-IRl  A delta  pitching  moment  coefficient  to  be  added  for 
IBLADE  = 3. 
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Output : 


CL  Blade  section  coefficient  of  lift 

CD  Blade  section  coefficient  of  drag 

CM  Blade  section  coefficient  of  moment 

1.3.3  BIRDl 

The  subroutine  BIRDl  has  the  entry  points  BIRDl,  PBIRD,  ABIRD,  BBIRD  and 
COR.  The  code  by  in  large  computes  and  integrates  the  isolated  control 
gyro  equations  of  Volume  I,  Section  6.11. 

Entry  BIRDl  computes  the  isolated  gyro  rotating  and  ncnrotating  inertias. 
The  ZDD,  ZD,  Z array  is  zeroed.  The  array  grouping  contains  the  shaft 
bending  and  gyro  degrees  of  fx'eedom.  Shaft  bending  accelerations  are  com- 
puted in  ACCEL  and  integrated  within  the  ZDD  array.  The  gyro  equations 
are  developed  totally  within  the  ZDD  context.  A summary  of  the  array  is 
given  in  Table  1-5. 


The  entry  PBIRD  computes  the  isolated  control  gyro  cyclic  actuators  time 
» constants.  This  entry  also  equivalences  the  shaft  bending  accelerations  to 

- ZDD  format.  Adams-Bashforth,  open  form  integration  for  the  ZDD,  ZD,  Z sets 

k completes  PBIRD. 
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The  control  flap  feedback  of  Volume  I,  Section  6.11.1^  is  coded  in  entry 
ABIRD.  This  code  is  directly  identifiable  with  the  symbols  in  the  text 
except  for  GLFEED,  GMFEED,  the  roll  and  pitch  feedback  moment  components 
respectively. 

The  final  isolated  gyro  equai-lons  of  Volume  I,  Section  6.11.3  are  assem- 
bled in  the  entry  BBIRD.  This  coding  can  be  traced  from  FKDD  = ♦q  and 
FDDD  = 'e  . Entry  COR  does  the  bookkeeping  for  the  values  of  old  accelera- 
tions and  velocities  required  by  the  integration  scheme.  The  term  NZ  is 
the  ZDD, etc.,  variable  size.  It  is  2 for  cases  without  shaft  bending  and 
U with. 

1.3. Ji  BMOVE 

BMOVE  is  a bookkeeping  subroutine.  The  prime  function  is  to  update  the 
blade  bending  modes  V(J,  (l,  2,  3).>  I)  from  the  YDD,  YD,  Y arrays.  The 
index  J is  the  mode  number  (l  to  3).  The  inner  index  (l,  2,  3)  is  for 
position,  velocity,  accele-^ation  respectively.  The  last  index,  I,  is  the 
blade  number.  The  array  FDD,  PD,  P corresponds  to  YDD,  YD,  Y for  the 
fourth  mode  per  blade.  The  index  of  the  "P”  arrays  is  1 to  which  is 
equated  to  1*,  8,  12  and  l6  of  the  "Y"  array  terms. 

A special  version  of  the  "P"  arrays  uses  only  information  from  ''Y"  array 
from  blade  i*.  This  data  is  distributed  to  the  "P"  locations  to  form  a 
reactionless  pitch  horn  bending  mode. 

1.3.5  BSCALE 

This  subroutine  figures  an  appropriate  plot  scede  factor  for  the  plotted 
variables  to  give  values  of  2,  5 or  10  per  line  of  plot  output. 

1.3.6  CMLOOK 

CMLOOK  is  the  executive  subroutine  for  the  seven  table  lookup  scheme  moment 
data  acquisition.  This  subroutine  calls  the  interpolation  routine  XTRP 
which  operates  on  the  moment  data  sets  T2CM,  TISCM  and  TICM.  The  moment 
data  selection  is  developed  in  Section  i».l  below  and  Volume  III,  Section 

3.3.7. 3. 


1.3.7  CNTROL 

The  subroutine  CI7TR0L  sets  the  command  control  positions  in  FLY.  It  con- 
sists of  the  entries  CNTROL,  APCNTL,  PICNTL.  Entry  CNTROL  equates  input 
command  time  steps  to  RA  inputs  and  initializes  constants  used  in  the  sub- 
routine. The  entry  APCNTL  in  conjunction  with  the  subroutine  PILOTA  for  a 
profile  follower  autopilot.  Insufficient  information  is  available  to  docu- 
ment this  section,  and  it  is  currently  not  considered  operational.  Command 
cyclic  stick,  main  rotor  and  tail  rotor  collectives  are  computed  in  entry 
PICNTL.  The  input  se..  of  position  vs.  time  (20  points)  is  interpolated  for 
the  computation  time  point.  A set  ~f  stability  augmentation  systems  are 
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meshed  with  these  calculations.  These  devices  are  noted  in  the  code  and 
diagrammed  in  Volume  III,  Section  3. 3. ^.2. 

1.3.8  COULOM 


Coulomb  friction  is  calculated  by  the  subroutine  COULOM.  COULOM  is  used 
for  swashplate  and  feather  bearing  friction  as  explained  in  Volume  I, 
Section  6.10  and  Fig'jre  6-3. 

1.3.9  DERIV 

This  subroutine  codes  the  equations  of  Volume  I,  Section  5.5«1  at  entry 
DERIV  and  the  equations  of  Volume  I,  Section  5. 5*2  at  entry  DERIVF.  In 
DERIV  the  earth  Euler  angles  are  IHETE  (6  ) and  PHIE  The  hub  axis 

accelerations  without  gravity  are  UHD,  ’/HD,  WHD.  These  terms  correspond 
to: 


i 


V 


(1-1 r) 


The  earth  angle  Euler  rates  sire  computed  as  PHIDE,  THETDE,  PSIDE. 

Entry  DERIVF  computes  the  fuselage  acceleration  without  gravity  (UPD,  VFD, 
WFD)  and  the  fuselage  velocity  (UF,  VF,  WF).  Txie  gravity  vector  in  fuse- 
lage coordinates  is  GFX,  GFY,  GFZ.  Ihe  shaft  bending  angles  ijfg,  0g  are 
coded  as  PHIS,  THTS. 

1.3.10  DWASH 


This  subroutine  is  a set  of  numerical  first-order  time  lags  applied  pri- 
marily to  main  rotor  and  tail  rotor  downwash  terms.  Derivatives  are  also 
calculated  from  time  step  increments  in  the  doxTiwash  components  and  the 
respective  time  constants.  The  subject  is  developed  in  Volume  I,  Section 
T.2.2.1*.  The  main  rotor  downwash  terms  are  wTMR,  PIMR,  QIMR.  These  ele- 
ments correspond  to  pi^p,  qi^p*  The  derivative  counterparts  are 

WIMRD,  PIMRD,  QIMRD.  The  tail  rotor  downwash  is  WITR. 

The  tail  rotor  flapping  dynamics  are  handled  in  the  same  manner  as  the 
quasi-static  pitch  horn  bending  in  Volume  I,  Section  6.6.6.  The  period 
of  oscillation  of  the  true  dynamics  is  replaced  by  a first-order  lag  with 
a time  constant  equal  to  the  period.  The  tail  rotor  flapping  is  coded  as 
AITR. 
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1.3.11  ETORQl 


ETORQl  computes  the  engine  model  equations  of  Volume  I,  Section  6.12.2. 

The  first  entry  ETORQl  initializes  the  gas  generator  speed  storage,  NGPRM, 
NG.  The  maxiiaum  available  torque,  MZZEDX,  is  calculated  from  the  maximum 
horsepower,  ENGHPX  and  nominal  rotation  speed,  0. 

The  equations  are  coded  in  entry  ETORQ.  is  coded  by  noting  that  the 

acceleration  is  the  difference  of  velocities  oetween  time  points  divided  by 
the  time  increment  DT.  The  velocity  here  is  the  rotational  speed. 


4' 

ENG 


= NGPRM 


The  engine  perturbation  torque  equation  is  coded  as  follows: 


ENDMZZ 

MZZTRM 

PQENG 

PQEOM 


M. 


^NG 

'\m. 


TRIM 


Limiters  are  provided  to  prevent  the  engine  supplied  torque,  ENDMZZ,  f : '"n 
being  negative,  or  exceeding  the  maximum  available  torque,  MZZEDX. 

1.3.12  m 

FLY  is  primarily  an  executive  and  outputting  subroutine.  MAIN  calls  FLY 
when  TRIM  is  complete  and  this  subroutine  controls  the  program  until  the 
maneuver  is  complete,  at  which  time  control  is  handed  back  to  MAIN  for  the 
next  case.  The  calling  sequence  in  FLY  is  given  in  Figure  1-2.  The  dia- 
gram gives  a brief  description  of  the  major  functions  of  FLY  and  the  sub- 
routine it  calls.  Each  entry  into  each  subroutine  is  enclosed  by  a 
separate  box.  Each  box  is  headed  by  the  entry  name  with  the  subroutine 
name  following  in  parenthesis. 

FLY  starts  by  initializing  itself  and  calling  the  initialization  entries 
CONTROL  and  PILOTI.  For  the  first  time  point  only,  the  program  then  Jumps 
over  a section  related  to  advaiwing  time  for  computation  which  need  not  be 
performed  at  zero  maneuver  time.  This  " otion  includes  calls  to  PRED, 
IGYRO,  DERIV,  DERIVF,  BMOVE,  PHORN,  TO  .,  ETORQ,  DWASH,  APCNTL  and  PICNTL. 
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The  reactionless  inplane  excitation  twang,  the  advancing  of  time  and  azi- 
muth, the  updating  of  spatial  angles  and  geometry,  and  the  updating  of 
downwash  save  veiriables  are  also  bypassed  for  the  first  time  frame. 

Next,  the  program  proceeds  to  update  the  collective  and  cyclic  main  rotor 
angles  as  function  of  the  svashplate  motions.  The  swashplate  motions,  in 
turn,  are  determined  either  directly  from  the  stick  or  from  the  swashplate 
control  actuators.  Then  another  series  of  subroutines  are  called;  LOAES, 
ACCEL,  DERIV,  CORECT,  BMOVE,  ABIRD  and  COR.  Next,  an  exit  test  is  made  for 
excessive  rotor  loads  or  the  end  of  the  maneuver.  All  the  rest  of  the  pro- 
gram is  devoted  to  collecting  plot  data  on  an  external  storage  device  and 
in  tabulating  output  data.  At  almost  the  end  of  the  subroutine,  there  is 
a section  for  initialization,  if  time  histories  to  an  expanded  time  scale 
are  desired,  for  an  extra  half  second  following  the  normal  time  histories. 
Every  time  point  is  plotted  for  the  expanded  plots,  whereas  normally  a 
number  of  tine  points  are  skipped  for  plotting  in  the  interest  of  keeping 
the  plot  length  reasonable. 

The  definitions  of  input  qmntities  equivalenced  to  the  input  RA  can  be 
found  in  Volume  III,  Section  3.2.  The  definitions  of  plot  parameters  are 
given  in  the  beginning  of  subroutine  RCPLOT, and  the  definitions  of  tabu- 
lated values  may  be  found  in  Volume  III,  Section  3. 3. 5 *2.  The  definitions 
of  key  variables  from  the  remaining  FORTRAN  variables  follow: 

AM  - number  of  time  points  during  normal  or  expanded  time  liistories 

AZIMUTH  - number  of  time  points  in  FLY 

TIME  - time  in  FLY 

PTIME  - plot  time 

PLATA  - time  history  parameters 

JSTEP  - time  point  counter  for  ZDD  and  ZD  integration 

KSTEP  - time  point  counter  for  YDD  and  YD  integration  for  the  blade 
degrees  of  freedom 

LSTEP  - time  point  counter  for  integration  of  the  remaining  YDD  and  YD 
variables 

IB  - number  of  blades  (U  presently) 

NVAh  - size  of  generalized  mass  matrix 

SCALE,  TSCLE  - revolutions  and  seconds  per  inch  of  plot  paper 
NSAVE  - number  of  time  points  to  next  plot  point 


NPTS  - number  of  points  plotted 

IPTB  - time  point  counter  for  plot  points 

NFREQ  - flags  expanded  time  histories. 

1.3.13  HARMl 

The  subroutine  HARMl  harmonically  analyses  selected  parameters  during  extra 
revolutions  added  to  the  TRIM  process.  Seven  harmonics  are  analyzed,  the 
first  being  the  "OP"  or  mean  value.  The  formulation  is  given  in  Vol- 
ume III,  Section  3.3.6.  HARMl  initializes  the  sine  and  cosine  for  the 
incremental  azimuth  angle  in  each  harmonic. 

The  entry  HARM  evaluates  the  integrals  required  in  the  coefficient  defini- 
tions. This  entry  is  used  every  integration  time  step  during  the  last 
TRIM  revolutions.  At  entry  HARM2,  the  fined  coefficients  eire  assembled, 
the  vector  form  of  the  analysis  is  computed,  and  the  harmonic  analysis 
printout  is  performed.  HARM2  also  adds  labels  to  the  tabulation. 

PESO  and  FBCO  compile  the  sine  and  cosine  component  of  each  parameter. 

PHA  and  AMPL  are  the  vector  phase  and  magnitude  representations.  SB  and 
CB  denote  the  sine  and  cosine  of  the  harmonic  signal.  HACYC  is  the  number 
of  rotor  revolutions  to  be  analyzed,  usually  one. 

1.3.14  IG^ 

The  purpose  of  subroutine  IGYRO  is  to  provide  the  capability  to  subdivide 
the  basic  integration  step  size  when  applied  to  the  gyro  degrees  of  freedom. 
The  relationship  of  interest  is: 


■*‘0YH0  ' 


(1-18) 


Where  GINT  is  an  input  quantity  RA(6l)  and  GINT  2:1. 

Note  that  upon  return  to  subroutine  FLY,  the  gyro  variables  will  have  been 
integrated  to  values  corresponding  to  a time  of  t + dt. 

Although  the  gyro  and  flexible  shaft  degrees  of  freedom  are  integrated  in 
the  common  routine  BIRDl,  the  flexible  shaft  and  subinterval  integration 
step  options  are  not  compatible.  This  is  because  the  gyro  degrees  of  free- 
dom are  explicitly  solved  equations,  and  therefore,  do  not  depend  on  the 
mass  matrix  operations.  Thus  the  degrees  of  freedom  are  not  affected  by 
the  step  size  change.  Such  is  not  the  case  with  the  flexible  shaft  degrees 
of  freedom  which  operate  through  the  normal  channels.  Thus  they  are 
directly  tied  to  the  mass  matrix  calculation  timing. 
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1.3.15  INTO 


Subroutine  INTO  is  the  system  integration  routine.  The  integration 
algorithm  used  is  the  four-point  Adams-Bashforth  open  formula 


The  subscript  n identifies  the  present  time  poin-*  , and  n+1,  n-1,  n-2  and 
n~3  neighboring  time  points  newer  and  older  than  n.  dt  is  the  increment 
between  time  points. 

Multipoint  formulas  such  as  this  require  companion  starter  formulas  to 
accumulate  back  values.  The  one  used  in  REXCR  is  the  simple  Euler  formula 


This  subprogram  has  three  entry  points.  INTO  is  the  routine  initialization 
entry.  Old  values  of  the  velocities  and  accelerations,  YDl,  YD2,  YD3, 

YDDl,  YDD2,  YDD3  are  zeroed  here.  Entry  PRED  evaluates  the  appropriate 
integration  formulas  and  advances  time. 

Entry  CORECT  performs  bookkeeping  functions.  Back  value  arrays  are  shifted 
to  prepare  for  another  integration  step.  Note  that  programming  is  pro- 
vided to  make  the  first  time  point  in  FLY  be  the  same  as  the  last  point  in 
TRIM.  This  is  done  by  back  stepping  the  integration  one  inteival  at  the 
end  of  TRIM  so  that  the  first  operation  in  FLY  will  bring  the  integration 
up  to  date. 

Since  the  integrators  are  not  restarted  fron  the  changeover  from  TRIM  to 
FLY  an  anomaly  can  occur.  If  the  integration  step  sizes  for  the  two  modes 
are  different,  then  the  results  for  the  first  few  points  into  FLY  are  not 
strictly  correct  until  the  old  integrators  are  filled  with  FLY  generated 
values . 

1.3.16  LOADSl 


The  entry  point  LOADSl  p<!rforms  a number  of  initialization  and  bookkeeping 
operations.  The  fuselag-s  lift,  drag,  moment  table  argument  (ALFA)  is  con- 
verted to  radians  and  a number  of  storage  locations  are  zeroed  out.  Also 
tail  rotor  aerodynamic  coefficients  are  precomputed. 

Main  rotor  coefficients  jre  precomputed  as  the  first  item  of  entry  LC'DS. 

An  array  of  main  rotor  aerodynamic  loads  in  rotating  coordinates  are  com- 
puted in  the  vector  FR.  FR  is  built  frcm  the  F array  from  sweep.  The  sec- 
ond index  is  blade  number.  Examples  of  the  ter.ms  are  vertical  load,  FR(9), 
roll  moment  referenced  to  blade  Hi,  FR(lO),  and  pitch  moment  referenced  to 
blade  #1,  FR(ll). 
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The  rotating  components  are  resolved  to  nonrotating  components  (MR)  by  the 
Euler  rotation  4'p.  The  terms  SCY  and  CCY  are  the  sine  and  cosine  values 
of  4'p.  These  loads  are  used  in  ACCEL  to  build  up  the  generalized  forces 
as  well  as  main  rotor  induced  velocity  calculations. 

The  new  values  of  Wj_jj,jpj  (WIMRN),  PiWR  (PIMRN)  and  q.it/TR  (QIMRN)  aro  computed 
as  given  in  Volume  I,  Section  J.Z.2.2.  Next  the  tail  rotor  loads  are 
partially  calculated  as  given  in  Volume  I,  Section  7.5.1. 

Wake  angle  logic  determines  if  rearward  flight  or  hovering  conditions  exist 
to  avoid  obtaining  an  improper  answer  from  the  use  of  the  iu*ctangent  func- 
tion. Using  the  evaluated  wake  angle  (Al),  the  wake  angle  tables  are 
interpolated  for  the  interference  factor.  Note  the  FXTN  interference  table 
contains  both  fuselage  and  tail  location  factors.  For  1=1  the  table  is 
used  to  get  the  fuselage  factor  FX.  For  1=2,  the  table  produces  the  tail 
wake  factor  which  is  used  in  the  temporary  computation  registe.'*  THBODY. 
Vertical  velocity  terms  WBODY  and  WBODYD  are  then  computed  using  this 
interference  data. 

•An  array  of  derivative  fuselage  (non  main  rotor)  aerodynamic  loads  are 
computed  as  FNW  and  are  directly  identified  with  Voliune  I,  Section  7.**.1. 
Static  data  are  added  to  these  terms.  The  0^,  Cp,  and  are  interpolated 
as  a function  of  local  angle  of  attack  ALFAB  and  input  data  ALFA(I).  The 
combined  FNW  array  is  resolved  to  fuselage  axis  as  per  Volume  I,  Section 
7.J*.1  in  array  FN.  The  tail  rotor  loads,  FTR,  and  propeller  table  lookup, 
FP,  follow  Volume  I,  Sections  7.5.2  and  7.6.1.  The  final  operation  is  to 
sum  FN,  Fl'R,  and  FP,  and  convert  to  the  actual  air  density.  This  final 
.summation  is  stored  as  QLOADS  for  use  in  ACCEL. 

1.3.17  LSTDAT 


This-  subprogram  is  called  once  at  the  outset  of  a computer  run  and  per- 
foms  two  major  functions.  One,  it  provides  a card  image  listing  of  the 
data  deck  Just  as  it  was  submitted.  This  provides  an  exact  record  of  the 
data  for  the  given  run.  Second,  LSTDAT  prepares  a working  data  set  by 
transferring  the  data  to  an  I/O  unit  NUNIT.  NUNIT  is  currently  assigned 
the  number  3.  See  Volume  III,  Section  6. 

1.3.18  MAIN 

As  the  neime  implies,  th5,s  x’outine  controls  the  REXOR  computation  sequence. 
The  firsc  function  is  to  initialize  the  problem  through  internal  opera- 
tions and  calling  of  the  initializing  entries  of  the  subroutines  LOADSl, 
SWEEPl,  ACCELl,  INTO,  TORSI,  BIRDl,  BMOVE,  ETORQl.  Control  then  passes  to 
TRIM.which  in  fact  is  an  executive  routine.  On  completion  of  a successful 
trim  or  time  limit,  control  transfers  to  FLY.  If  a number  of  input  cases 
are  being  run,  control  cycles  back  to  label  10  for  additional  runs. 
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1.3.19  MIC  26 


Within  REXOR,  the  equations  of  motion  are  stated 

|qt  = jAq}  + jq^g^J  (l-2l) 


where 

)A*q(  = [m]"^  |af| 


(1-22) 


M and  A?’  are  given.  Aq  can  be  computed  by  first  inverting  M,then  per- 
forming the  indicated  multiplication.  A more  efficient  method  is  to  solve 
for  the  product  directly  by  solving  the  linear  system 


[m]  |Aq(  = {af}  for  jAq( 


(1-23) 


It  is  further  known  that  the  mass  matrix,  [m],  is  positive  definite  and 
symmetric.  Subprogram  MIC26  is  a general  algorithm  for  the  solution  of 
simultaneous  equations  of  the  form 


[A]  {xj  = 


n-sh) 


where  the  coefficient  matrix  is  positive  definite,  symmetric.  The  algo- 
rithm is  a Cholesky  decomposition  of  [a] , followed  by  a forward-backward 
substitution.  The  algorithm  is  presented  below. 

Cholesky  Method  for  Symmetric,  Positive  Definite  Matrices . 

Theorem:  Let  A be  symmetric,  positive  definite.  Then  A can  be  factored 

in  the  form 


T 

LL  = A 


(1-25) 


where  L is  a lower  triangular  matrix  (i.e. , L = (t.,)  where  £.  = 0 for 

j > i) 
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Cholesky  Method:  Let  A be  n x n,  symmetric,  positive  definite 


A = (ajj) 


a,  . = a,, 
ij 


(1-26) 


Assume  A is  factorable 


Then, 


II 


Then 


A = LU 


(1-27) 


a . . = > i , u,  . 

ij  ij  kj 

k=l 


(1-28) 


U = L" 


’ \j  ^jk 


(1-29) 


n 

“ij  " I Hk  »Jk  i = 1.  ni  J = 1,  n (1-30) 

k=l 


But  L is  lower  triangular,  which  implies  = 0 for  k > i.  Therefore, 


i 

'IJ  ° 2 «ik  »Jk  i = 1.  n;  3 = i.  n (1-31) 

k=l 

Equation  (l-3l)  forms  the  basis  of  the  decomposition.  The  elements  of  L 
are  found  as  follows. 
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^11  = / ^11 


(1-32) 


^IJ  ""  ^11  ^Jl 


^jl  ~ ^lj/^11  * j = 2,  n 


(1-33) 


“u  = S »ik"  = S ii/  * ^ii' 

k=l  k=l 


for  i = 2,  n can  be  solved  for  il . . . 

li 


i .,  = 
li 


= V “u-I 


U-3i*) 


‘ij  “ I hi 


LJ  Z-  ^ik  'Jk  *'ii  ''ji 
k=l 


"'Ik  •"  ^n  > i = 2,  n (1-35) 


^ik  ^jk)  /hi 


(1-36) 
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With  the  decomposition  of  A,  the  system  (l-2l<)  can  be  solved  as  follows. 
Substitute  A = LL^  into  (2).  LL^x  = b 


- tT 
g = L X 


Solve  Lg  = b by  forward  substitution.  Namely, 


'i  “ - 2 »ik  %)/‘ 


i = 2,  n 


Finally,  find  x from  (l-37)  bj  backward  substitution. 


X = g /t 
n ®n  nn 


n 

'‘j  " (®J  ■ S ^kj  O/^jJ 


k=j+l 


Equations  (l-32)  through  (1-36)  and  (1-38)  through  (1-^*1)  form  the 
algorithm. 

The  arguments  in  the  subprogram  cailling  sequence  are 


MIC26  (A,  B,  N,  M,  L) 


whe  e A is  the  coefficient  matrix  of  a maximum  dimension  M x M.  A sub- 
matrix problem  of  size  N x N can  be  solved.  B is  the  right-hand  matrix  of 
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dimension  N x L on  entry  and  the  solution  matrix  on  exit.  The  algorithm 
assumes  A is  symmetric  and  only  operates  on  the  lover  triangular  portion, 
A is  destroyed  on  exit. 

An  error  exit  can  occur  from  the  routine.  While  forming  the  elements  as 
given  by  equation  (l-35)j  the  argvunent  of  the  square  root  function  is 
tested.  If 


I 


\ 

\ 

i 

i 


then  the  procedure  is  terminated.  This  is  interpreted  as  an  indication  of 
failure  to  be  positive-definite.  A normal  exit  from  the  routine  is  indi- 
cated by: 


A^^  = 0 (1-1*3) 


An  error  exit  is  indicated  by: 


Ajj  = S U-W) 


Where  j is  the  row  number  at  failure. 

The  subprogram  MIC26  is  computed  in  double  precision  on  IBM  hardware.  This 
means  the  matrices  A and  B are  double  precision.  The  CDC  version  is  single 
precision.  See  Volume  III,  Section  6. 

1,3.20  MINVR 

This  subprogram  will  compute  the  inverse  of  a matrix  by  Gaussean  Elimina- 
tion with  pivoting.  The  inversion  is  performed  in  single  precision  and 
the  matrix  dimensions  have  been  specialized  to  6 by  6.  The  calling  se- 
quence arguments  are  MIhTR  (A,  N,  B,  lER,  W) 

where 

A - input  matrix  of  order  N x N N - 6. 

B - output  matrix  A ^ 
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W - work  array  of  size  6 by  6 
lER  = error  indicator 

0 - normal 

1 - pivot  element  zero,  (no  inversion) 

MINVR  is  called  only  once  from  TRIM  and  only  if  one  of  the  unsupported 
trim  options  is  activated.  Its  sole  purpose  is  to  provide  a constant 
sensitivity  matrix  for  a Nexrton-Raphson  type  of  iterative  equation.  Given 
A,  find  B = A~1  to  be  used  in  an  equation  of  the  form 

Vl  = \ ® 


where  x is  a control  vi. ctor  and  f is  a trim  error  vector. 
1.3.21  MPRNT 


MPRNT  is  a generalized  matrix  formatted  print  routine.  It  is  called  by 
subprogram  ACCEL  to  provide  a print  of  the  mass  matrix.  The  printing  is 
optional.  If  activai,ed,  via  input,  the  mass  matrix  is  printed  once  at  the 
beginning  of  TRIM  and  FLY. 

1.3.22  PLATE 

Subprogram  is  a dummy  date  routine  which  is  to  be  supplied  by  the  installa- 
tion of  interest.  The  calling  routine  displays  the  date  to  be  presented 
as  a 1 word,  8 character  literal. 


CALL  PLATE  (NLATE) 


The  argument  NLATE  is  typed  as  REAL  *3  for  lEM  software.  The  image  is  ex- 
pected to  be  of  the  form 


MO  - LA  - YR 


The  CLC  counterpart  is  typed  real  and  is  also  of  the  form. 


MO  - LA  - YR 


giving  eight  characters  left  adjusted. 
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i 1.3. 23  PHORN 

This  subroutine  solves  the  isolated  pitch  horn  bending  er mtion  of  Vol- 
ume I,  Section  6.6.6.  TPH  is  Xp  . The  term  TFA  corresponds  to  blade 
feathering  moment  Mp^.  DPF  represents  and  DPFD  gives  ^p^  p^j. 

1.3.21+  PRINTl 

I Subprogram  PRINTl  has  multiple  entry  points  all  called  from  subprogram 

I READIN.  Entry  PRINTl  presents  a formatted  printing  of  the  master  data 

deck.  Each  input  is  identified  by  listing  its  address,  program  name, 

I ■ and  value.  Tabular  information  is  presented  as  such  with  some  grouping 

I for  easy  identification.  Entry  PRINT  gives  case  data  as  an  exception 

\ report.  Case  data  which  is  different  in  value  from  the  master  data  value 

' is  presented. 

1.3.25  PROP 

PROP  is  a data  bank  of  propeller  tables  for  the  AH-56A.  This  data  is  used 
by  the  subroutine /entry  LOADSl/LOADS  via  the  bivariant  interpolation 
ZTRP2.  The  first  grouping  of  data  PROPl  contains  C,j,  information  as  a func- 
tion of  blade  angle,  6,  and  advance  ratio,  J.  The  second  grouping,  PR0P2 
contains  similar  information  for  Cp. 

1.3.26  RCPLOT 

t 

This  subprogram  performs  all  the  CALCOMP  plotting.  The  coding  has  been 
arranged  to  give  optimal  results  when  used  with  the  12-inch  drum  CALCOMP 
machines.  Other  units  may  require  some  reprogramming  to  achieve  useful 
results.  The  routine  performs  as  described  in  Volume  III,  Section  3.3.5. 

It  is  called  optionally  from  TRIM,  for  trim  and  harmonic  analysis  plots, 
and  from  FLY  for  maneuver  time  histories.  It  will  also  optionally  prepare 
a data  set  of  signals  for  transform  analysis.  The  calling  sequence  is  as 
follows : 


RCPLOT  (NVAR,  NVEC,  NPTS,  CASE,  PLTINC) 


where: 

NVAR  - number  of  parameters  to  be  plotted. 
NVEC  - parameter  code  table. 

NPTS  - number  of  points  to  plot. 
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CASE  - case  number. 


PLTINC  - time  between  data  points. 

NVEC  = NVECl  = RA(301)  through  RA(3l*0)  in  TRIM,  or  NVEC  = NVEC2  = RA(i801) 
through  RA(1060)  in  FLY.  Up  to  90  parameters  can  be  plotted,  however  num- 
bers greater  than  6o  fold  over  on  a modulus  60  basis. 

Some  input  factors  and  controls  are  transmitted  through  the  COMMON  state- 
ments. The  scaling  is  done  by  the  subprogram  BSCALE,  if  automatic  scaling 
is  used.  SVEC  = RA(l85l)  through  RA(l900)  sets  the  scale  factors  in  FLY. 

1.3.27  READIN 

This  subprogram  processes  the  input  data,  building  cases  as  defined  by  the 
control  cards.  The  master  data  configuration  is  preserved  in  array  RAS. 
The  working  case  configuration  is  defined  in  array  RA.  READIN  senses  the 
end  of  a computer  run,  terminates  the  plot  data  file  and  exits  execution. 

1.3.28  STALL 


STALL  computes  the  elements  of  dynamic  stall  given  in  Volume  I,  Section 
7.2.3.**.  The  subroutine  is  called  from  SWEEP,  and  is  operative  with  the 
fast  aero  lookup  scheme  (Volume  I,  Section  7 .2. **.3  and  Section  *».2  of  this 
volume).  The  partial  derivative  (3Cl/3«)  is  calculated  by  two  ceills  to 
AER'i  with  angles  of  attack  0 and  0.01  radian.  The  CL  differences  are 
multiplied  by  100  and  termed  DCLDAO.  The  dynamic  stall  delay,  y*  is  com- 
puted as  TEMP  with  Mach  limitations  of  0.25  and  0.6.  The  reference  angle 
of  attack  Opgp  is  then  calculated  as  AD. 

A logic  sequence  given  in  Volume  I,  Section  7.2.3.**.l  is  used  to  select 
the  proper  (8Cj^/3a)  so  that  the  static  lift  slope  line  is  not  exceeded. 

The  value  selected  then  occupies  TEMP.  CLR  is  the  completed  lift  slope 
value. 

The  dynamic  stall  moment  calculations  are  produced  simultaneously;-  with  the 

above  operations,  and  in  accord  with  Volume  I,  Section  7.2. 3. **.2.  The 

factor  K is  FACTM  and  is  ADM.  The  final  moment  coefficient  is  CMR. 

nhr 

1.3.29  SWEEPl 


SWEEP  is  the  key  routine  to  the  REXOR  program.  Elements  of  Volume  I,  Sec- 
tions 5.5.5  and  6.6  are  computed  here.  The  blade  equations  of  motion 
ingredients  are  computed  for  an  azim-ath  step  for  each  blade  and  radial 
station. 

SWEEPl  has  the  entries  SWEEPl  and  SWEEP.  A large  amount  of  information 
is  equivalenced  from  input  RA  set  in  SWEEPl  as  well  as  out  to  names  used 
in  other  routines.  The  RA  names  can  be  matched  against  the  input  set 
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f 


> 


t 


I 


•J 


description  of  Volume  III,  Section  3.2.  If  the  fast  aero  option  is  speci- 
fied the  five-point  input  thickness  and  design  lift  coefficient  points  are 
interpolated  to  the  blade  radial  stations  used.  The  thickness  values  are 
TCN(K)  and  the  lift  coefficients  are  CLN(K). 

Radial  stations  are  set  for  output  functions  and  for  use  in  other  subroutines. 
These  are  NTHl,  NTH2,  NSTAF.  The  first  two  are  stations  to  read  out 
dynamic  torsion.  The  third  is  the  point  to  read  out  effective  sweep  and 
droop.  Program  logic  outputs  data  at  the  closest  blade  station.  Blade 
mass,  BLMASS,  and  rotor  polar  inertia,  IZZR,  are  precomputed  at  this  point. 
The  factor  DSS2(K)  is  the  coefficient  string  for  trapezoidal  integration. 


Coordinate  transformation  parts  and  feather  axis  data  are  next  precomputed. 
BO  is  B , GAM  is  Y , BF  is  B„,»  PHIREF  is  <)>„„.  Transformations  calculated 

O r A K&r 

are 


M = Kf  P v]' 


M = [\f  bJ 


(1-U6) 


(I-U7) 


and  partially 


(1-1*8) 


The  elements  XRI,  ZRI,  XRC,  ZRO  are  the  inboard  and  outboard  tension-torsion 
attach  point  locations. 

The  entry  SWEEP  is  called  from  LOADS  for  each  azimuth  step.  First  the  hub 

set  airflow  angles  are  computed  per  Figure  5-l<  of  Volume  I.  A2  is  02  and 

sine,  cosine  components  of  'jiy  are  SPSW  and  CPSW.  Main  rotor  ground  effect 

coefficient  f..,„  is  FHl. 
iMR 

Inboard  and  outboard  feather  bearing  displacements,  velocities,  and  accel- 
erations are  then  determined.  The  modal  coefficients  used  are  pre- 
interpolated to  the  bearing  locations. 
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These  coefficients  are  defined  as: 

FBLII,  first  inplane  made 
FBLIF,  first  flap  mode 
FBL2F,  second  flap  mode 
With  the  argument: 

(1  = inplane  deflection  1 = inboard  location  \ 

or  , ) 

2 = flapping  deflection  2 = outboard  location  / 


The  modal  variable  is  coded  eis: 


mode  number 


type 


I for 


V 


1,  2 or  3 


1 = displacement 

2 = velocity 

3 = acceleration 


blade 

number 


Using  the  data  calculated  above,  the  angles  Y'  (YPFA)  and  Z 


FA 

well  as  associated  time  derivatives  are  programmed.  Next, 

and  SPSI  set  up  the  cosine  and  sine  of  for  each  blade. 

a 


the^ 


(ZPFA)  as 
erms  CPSI 


The  remaining  elements  prior  to  integrating  (sweeping)  out  the  blade  set  up 
needed  transformations  and  partial  derivatives. 


Key  items  are: 


CPPW 

SPPW 


relative  wind  angle  resolved  into  BLn  set 


VKOROC  chordwise  velocity  with  ground  effect 


^1’  ^2’  ^3  acceleration  resolved  to  BLn  set 

XTl,...,  ZT3  angular  acceleration  cross  product  terms 


ANGR(l) 

ANGR(2) 

ANGRD(l) 

ANGRD(2) 


> BLn  set  roll  and  pitch  rate  and  acceleration 
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THF  true,  total  feather  command,  and  numerous  partial  derivatives: 
PFTHO  = 

PFAIS  = 

PFBIS  = 

PFHG(I)  = 3'J’p/3Zg 
PFGK(I)  = 

r u 

PFGD(I)  = 3^p/30g 
PFBLFB  = 

PFBLF  = 3<^bLe/'^F 
PFBLTB  = 3^bLe/^®H 
PFSI(I)  = a<^p/3>i<p 
PXRXB(I)  = pX^Ln/'x 
PYRXB(I)  = 3Y3j^/3x 

PZRZB(I)  = - 3^33^/ 

PYRYB(I)  = ^^BLn'^^Y 

PXRYBd)  = 

DELPHI  “ 4*_ 

F 

YPVl  = 3Y'  73A 
FA 

• • 

« • 

• • 

ZPV3  = 3Z’^^/dA^ 

TPn  = 3<^^/3A^ 


45 


Using  the  subroutine  XTERMO  and  multiple  entry  points,  these  transforma- 
tions are  calculated. 


H ■ W [vj  [v,j 

H ■ [xl  [’■■„]  [’■•  J 


(1-51) 


(1-52) 


(1-53) 
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FAJ 


il-5h) 


(1-55) 


(1-56) 


cT30],  [t36],  [t37],  [t39],  [T33]  are  partial  elements  of  double  dot  trans 
formation  r^j^  from  r^^g  data.  [t1j2]  is  the  transform  elements  less 
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translation.  The  preamble  concludes  by  zeroing  a number  of  storage 
locations . 


The  majority  of  the  calculations  within  SWEEP  eu*e  governed  by  two  FORTRAN 
DO  loop  indicies,  I and  K.  I is  the  blade  index  and  K is  the  station 
index.  Once  the  preamble  calculations  for  each  blade  are  determined,  the 
blade  spatial  integration  loop,  controlled  by  index  K,  begins. 


First, the  bending  acceleration  array  is  built  as 


veci6' 

VECU(l)' 

r -V 

0 

VEC17 

-■  = . 

VECl»(2) 

► = 4 

VEC18 

Ik 

VEC2*(3) 

BEND 


(1-57) 


The  modal  coefficients  are  coded  as  BMSII  for  first  inplane,  BMSIF  for 
first  flap,  and  BMS2F  for  second  flap.  These  vsiriables  use  the  following 
argument  scheme. 

1 = Y deflection 

2 - Z deflection 

3 = Y' slope 
k = Z* slope 

The  modal  variable  Vs  are  used  as  previously  described.  PHIT  is  the  in- 
stantaneous blade  twist,  and  is  used  to  complete  [Ti^S], 

Next,  the  transformation  of  Volume  I,  Section  5*5.5  from  BLE  to  BLn  axis  is 
developed. 

Waypoints  are 


K = blade  station. 


Y 

r > 

XSTAT 

BLE 

YSTAT 

* = * 

^BLE 

ZSTAT 

Z.  T ri 

..  * 

L JiLEj 

STATIC 


(1-58) 
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’ ' 
XBRA 

BLAD(l) 

YBRA 

BLAD(2) 

ZERA 

BLAD(3) 

foUMOd)! 

fo 

1 

f 

DUM0(2) 

DUMOO) 


J 

• > 

XD 

^LE 

YYD 

- = - 

^BLE 

ZD 

^^BLE^ 

BLE 


L ^BLE  J 


- < 


BEND 


LzJ 


INBOARD 


BLn 


X 

YY 


Li 

V.  J 


Y = A 


h 

Y 


leI 


BLE 


■"BLE  J 


BLn 


(1-59) 


(1-60) 


(l-6l) 


Likewise  for  accelerations  noting  the  X values  are  done  later. 


XDD 

’■^le' 

YYDD 

► = . 

^BLE 

ZDD 

.^ble^ 

BLn 


(1-62) 


Similar  procedures  are  followed  for  the  bending  slopes.  Again  noting 
waypoints: 


VECIO 

^ vEcu  y 

VEC12 


1 ^'ble 


BLE 


BEND 


(1-63) 
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! 


VEC13 

" % 
0 

VECT3(1) 

VECllt 

^'ble 

* = . 

VECT3(2) 

VEC15 

.^'ble. 

BEND 

.VECT3(3), 

VECT5(1) 

' -N 

0 

VECT5(2) 

► - , 

Y' 

BLE 

VECT5(3) 

?7  f 

BLE^ 

Going  through  coordinate  tiansforma'cions  gives 


- YP 

• = . 

^'ble 

.^’ble^ 

1 

r h 

BLn 


YPDI 
I^ZPDIj 


t < 


BliE 


BLE 


BLn 


I' 

’ ^ 

YPDDI 

« = , 

Y' 

BLE 

ZPDDI 

w w 

^ bleJ 

BLn 


(1-65) 


(1-66) 


(1-67) 


(1-68) 


A series  of  transforms  are  then  constructed.  These  elements  are  used  to 
calculate  angular  transforms.  First,  the  elements  of  TZB  are  rT„,  1'^ 
Other  pairs  are:  TZBI  with  TT-,  1,  TYB  with  rT„,  "I  [ BENDJ 

and  TYBI  with  [t  , 1 L BENDJ  L BENDj 

L BENDj 
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Using  these  formulations  leads  to 


> BUS 


(ANGBLD)  = < 4 


^^BLe/^^1,2,3  “ 

. 

r 

4 


YDPARl 


I YDPAR2 


YDPAR3 


(1-73) 
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I ZDPARl 
ZDPAR2 
ZDPAE3 


f 

► 

Jaz* 

1 BLE' 

1 

' > 

wU 

P 

- W5 

• = ■ 

q,  ' 

W6 

r 

V ) 

r 

W7 

p 

jw8 

1 

1 

► = ^ 

• 

q 

W9 

• 

r 

V 

BLn 


BLE 


(I-7I4) 


(1-75) 


(1-76) 


To  this  point  a number  of  Y and  Z components  have  been  developed.  Elements 
to  proceed  with  similar  X calculations  are  done  starting  after  label  320. 
Wembers  of  the  procedure  are: 


SYS  running  sian  of  spline  length  S 

YHAOBL  Y.„. 

ONA 

Yf.'ABLS  running  sum  of  Y^j^^  in  BLn 

YNABLR  increment  of  Y^„.  between  K and  K-1 

ONA 

ZNABLS  running  sum  of  Zj^^  in  BLn 

ZNABLR  increment  of  Z,,.  between  K and  K-1 

NA 

YNAS  running  sum  of  in  BLn 

YNAR  increment  of  Y„.  between  K and  K~1 

NA 

W value  of  X without  corrective  tern 


X W - X„.  = X,,,  in  BLn 

NA  NA 


Similar  operations  then  are  performed  to  get  X 
by  the  D and  DD  in  the  notation. 


NA 


See  Section  5.5.5.10 
of  Volume  I 


j 

• • 

and  X„,  These  are  noted 
NA 
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isrivative  calculations  prior  to  aeroiynaaic 


/l=Ai  /1=  $ 

PAAl2=A2|2=0 
\3  = \3  = f 


An  example  is  PAAl(2). 


PXBLPB 


PYBLFB 


PZBLFB 


PXBLTB 


PYBLTB 


PZBLTB 


PTBLSI 
PYBLSI  . = 
PZBLSI 


a^p  Vble 


“ ' «BLE/’*ln 


XV2  = 9, 


^LE*^^^2n 


"^3  = 
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At  this  juncture,  sufficient  information  exists  to  develop  the  aerodynamic 
loading  functions  and  complete  the  blade  calculations  of  this  subroutine. 
First  looking  at  air  velocities: 


UTl 

X 

UT2 

► = < 

Y air 

velocities 

UP 

Z 

BLn 


(1-81) 


'V2(D  ' 

X 

V2(2) 

V2(3) 

Y air 

velocities 

Z 

w * 

BLE 


Note  the  more  traditional  usage  is  equivalenced  at  this  point. 


(1-82) 


UC  - air  velocity  chordwise 

UN  - air  velocity  normal  to  chord. 

From  this  data  the  angle-of-attack  information  and  assorted  aero  rerms  are 
computed. 

ALFAR  - angle  of  attack 

/^LFDOT  - rate  of  angle  of  attack 

QC  - aerodynamic  pressure  times  chord 

RFREQ  - reduced  frequency 

AYAW  - spanwise  flow  angle  of  attack 

Depending  on  the  option  selection,  the  aerodynamic  loading  comes  frcm  call- 
ing STALL,  AERO  or  the  pair  XTRPU,  CMIX)OK.  The  first  candidate  is  used  for 
dynamic  stall.  Dynamic  stall  uses  fast  aero  (AERO)  which,  as  indicated, 
may  be  used  without  dynamic  stall  directly  from  SWEEP.  The  other  loading 
data  source  is  the  7 table  aero  lookup.  Here  the  subroutine  XTPRU  produces 

C..  and  C„.  C.,  data  comes  from  CMI/  K.  The  resulting  aero  loadings  are: 

L DM 


FNO 

normal  force 

FC 

. = > 

chordwise  force 

TX 

torsion  moment 

1 

BLE 


(1-83) 


The  ingredients  necessary  for  the  REXOR  equations  of  motion  are  assembled 
in  the  array  FI,  which  is  given  in  Table  1-6.  These  elements  are  per 
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TABLE  1-6.  FI  LIST  IN  SWEEP 


Element 

FI  (1) 


Use 

F^  3X/3A^  + F^  3Y/3A^  + F^  3Z/3A^  + 34>/3A^ 

+ 34'/8A^ 


FI  (2) 


See  F 

mn  FI  (3) 


Volume  I, 
section 
6.6. If 

FI  (If) 


FI  (5) 
FI  (6) 


F^  3X/3A2  + F^  3Y/3A2  + F^  3Z/3Ag 

+ 3(j)/  Ag  + 3'I'/3A2 

F^  3X/3A^  + F^  3Y/3A^  + F^ 

+ 3<I>/3A2  + 3<J^/3A2 


FI  (T) 
FI  (8) 
FI  (9) 


TABLE  1 

-6  - Continued 

Element 

Use 

FI(IO) 

^BLE 

Fl(ll) 

i 

BLE 

Aerodynamic  loading  only 

► 

FI(12) 

^BLE 

V.  J 

BL 

1 

FI(13) 

^BLE 

Fl(llt) 

My 

BLE 

Aerodynamic  loading  only 

► 

FI(15) 

L ^BLE  J 

BL 

FI(16) 

Aero  only  '' 

FI(IT) 

- 

, Diagnostics 

Fl(l8) 

BLE 

> 

FI(19) 

♦r*P 

Generalized  masses  completed  in 
sweep,  typical 

FI (20} 

A 

FI 

FI(21)  M. 
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TABLE  1-6  - Continued 

Element 

Use 

FI(22) 

"Vs 

FI(23) 

Vi 

Fl(2l*) 

^ (Typical) 

FI(25) 

M 

A A 

FI(26) 

Ma  a 
^2^2 

FI(2T) 

Ma  a 

Vs 

FI(28) 

^3^3 

FI(29) 

f, 

R 1 

FI(30) 

. 

Vi 

FI(31) 

^>1'  A 
R 3 

FI(32) 

^<1'  <li 
R^F 
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TABLE  1-6  - Continued 


Element 

Use 

FI(33) 

m aX/3A^ 

" V ^ 

Fl(3l*) 

m ax/aAg 

° \bI  '^2 

FI(35) 

m ax/aA^ 

" ^ A 

OBL  3 

FI(36) 

m 3Y/aA, 

X 

= ^ A 

OBL  1 

FI(37) 

m aY/aAg 

= ^ A 

OBL  2 

FI (38) 

m aY/3A^ 

“ “'OBI  S 

FI(39) 

m 3Y/3<l>p 

FI(1*0) 

m az/aA^ 

= M 

^OBL  ^ 

FI(1»1) 

m az/aAg 

= M 

7 A 

OBL  2 

FI{k2) 

m az/aA^ 

= M 

^OBL  S 

Fl(li3) 

m 3Z/3<f>„  = M 

^^OBL 

Fl(l*l4) 

”Vi 
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TABLE  1-6  - Continued 


Element  Use 


Fl(!l5) 

“Vs 

FI  (1*6) 

V3 

FI(1*7) 

“V. 

FI(1*8) 

“V. 

Fl(l49) 

Ve 

FI(50) 

“V3 

FI(51) 

“e  ♦ 

H F 

FT(52) 

(V.L 

FI(53) 

Fl(5l») 

m «{-X)  = M_  , 

'^OBL  PL 

FI(55) 

m*y=M^ 

OBL^BL 

FI(56) 

m »(-Z)  = 

^obl'^bl 

TABLE  1-6  - Continued 


Element 

Use 

FI(5T) 

" FI  (1)' 

FI(58) 

► Similar  to  < 

(2) 

but  aero  only 
(diagnostics) 

Fl(59) 

1 

Fl(60) 

Fl(6l) 

FI  (62) 

M,  , Uncoupled  dynamic 

V3 

FI(63) 

Mq  . Torsion  option 

Wp 

Fl(6h) 

^6  6 
PH  PH 

FI(65) 

F„  less  - Kg 

PH  PH*^PH 

term 

NOTES:  1.  FI's  summed  over  blade  and  stored  as:  F (I-65),  (n  •=  1 - k) 

2.  FI (66)  gives  remaining  part  of 

2 


Formulation  not  used  in  ACCEL. 
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radial  station.  The  array  F(J,  l)  stores  the  results  of  trapezoid  integra- 
tion of  FI.  The  index  J coincides  with  the  meaning  of  the  elements  of  FI. 

I is  the  blade  number.  The  F array  appears  in  LOADS  and  ACCEL. 

1.3.30  TRIM 


The  subroutine  TRIM  performs  the  function  of  producing  an  initial  set  of 
conditions  for  the  equations  of  motion  from  the  user  supplied  set  of  flight 
conditions.  TRIM  works  in  a number  of  different  modes  depending  on  the 
type  of  flight  condition,  such  as  free  flight  or  fixed  rotor,  and  the 
vehicle  configuration  itself.  Particular  items  are  the  arrangement  of  the 
lifting  surfaces,  and  the  control  system  configuration.  The  TRIM  options 
to  accommodate  the  different  required  modes  are  by  in  large  controlled  by 
CORAF  = RA(U2). 

TRIM  operates  as  a set  of  servo  loops  which  through  successive  iterations 
set  the  controls  such  that  the  desired  flight  path  is  followed.  The  pro- 
cediure  is  to  move  the  contr-- '.-motion  pair  in  question  to  null  the  differ- 
ence between  the  current  ca  'ulated  value  and  the  desired  final  value. 

Table  3-1  in  Volume  III  shows  the  control  or  variable  to  motion  pairing 
used  in  REXOR.  TRIM  does  not  compute  the  motion  of  the  entire  rotorcraft, 
but  only  computes  the  motion  of  the  rotor  assembly  feee  Table  1-T).  Other 
than  accelerations,  the  remt  ining  displacements  and  velocities  are  com- 
puted by  a set  of  static  relationships.  This  procedure  saves  time  by  con- 
straining the  flight  track  of  the  rotorcraft  during  TRIM,  and  eliminates 
what  would  otherwise  might  hi  a lengthy  chasing  procedure. 

TRIM  is  called  by  MAIN.  It  performs  executive  functions,  tabulates,  plots, 
specifies  motions  of  degrees  of  freedom  inoperative  in  trim, and  performs 
the  trim  algorithms  until  one  of  the  trim  tests  is  met  and  TRIM  exits. 

There  are  no  multiple  entries.  The  normal  trim  is  with  the  degrees  of 
freedom  operating  for  all  the  blades.  The  subroutine  includes  fast  trim 
procedures  using  one  blade  only.  These  trim  procedures  are  not  supported. 

The  program  initializes  the  new  case  down  to  label  9677.  At  this  point  the 
trim  loop  begins.  The  first  calculations  in  the  trim  loop  relate  to  updat- 
ing for  the  next  time  point  and  are  skipped  by  going  to  label  10^4  the  first 
time  into  the  trim  loop  at  time  zero.  Figure  1-1  illustrates  the  entire 
sequence  of  operations  in  TRIM  and  shows  the  initial  entry  from  MAIN. 

Beginning  with  label  104,  the  program  specifies  the  degrees  of  freedom 
inoperative  during  trim.  These  degrees  of  freedom  are  given  in  Table  1-7. 
The  equations  defining  these  degrees  of  freedom  are  taken  directly  from 
Volume  I,  or  are  derived  from  the  same  assumptions.  Note  that  in  TRIM  the 
fuselage  velocities  and  angles  are  the  prime  quantities  rather  than  the 
principal  axis  values.  This  scheme  is  used  to  make  the  trimming  quouLi- 
ties  closely  identifiable  with  normally  measured  quantities.  Hub  air 
velocities  may  be  expi'essed  in  terms  of  fuselage  data  by  inverting  the 
equations  of  Volume  I,  Section  5* 5*2; 
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TABLE  1-7.  SPECIFICATIOKS  FOR  DECREES  OF  FREIDOK  INOPERATIVE  DUHIHG  TRIK 
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The  quantities  (VT),  (ALPHA)  and  gp  (BET)  are  inputs  or  trim 
variables. 

The  remaining  item  is  the  pitch  attitude,  0g.  This  quantity  is  derived  as 
an  explicit  calculation  rather  than  the  looping  shown  in  Volume  I, 

Section  5-5.1  for  FLY. 

Starting  with  the  identity, 
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where 


[\] 


cos  a„  0 -sin  u^' 
n n 


sin  a„  0 cos  a„ 
n n* 


(1-87) 
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(1-91) 


= V /V 
h'  T 


(1-92) 


The  climb  angle  Y„  (GAMMA)  is  an  input  or  a trim  variable, 
n 

The  program  modifies  the  specifications  slightly  if  FLY  is  conducted  with 
the  shaft  fixed  simulating  a wind  tunnel  model  with  rigid  support  or  if  the 
pitch  and  roll  rates  are  specified  by  inputs.  Almost  all  the  specifica- 
tions are  complete  by  about  label  112, 

Once  the  trim  motions  are  specified,  the  prograir.  proceeds  to  call  a large 
number  of  subroutines  and  their  entries:  DERIV,  DERIVF,  LOADS,  ACCEL, 
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CORECT,  BMOVE,  and  ABIRD.  Note  that  if  harmonic  analysis  is  being  con- 
ducted, DERIV,  DERIVE  and  SWEEP  are  called  again  so  the  blade  loads  can  be 
recomputed  (aiid  made  timewise  correct)  using  the  accelerations  updated  in 
ICCEL  before  integrating  for  the  •"•elocities  and  displacements.  These 
calls,  along  with  some  fast  trim  programming,  occupy  the  program  to  just 
beyond  label  .10^+5. 

The  remaining  elements  partially  within  the  blade  radial  summation  loop 
(ends  at  label  399)  are  the  blade  loads  harmonic  analysis  (HSAVE)  and 
associated  terras. 

The  routine  ends  with  the  ^v,n°ion  torsion  pack  elements  aad  blade  feather- 
ing torque,  TIA.  The  blade  loop  ends  at  label  503. 

The  program  next  pr.-parts  output  for  harmonic  analysis  down  to  label  1053, 
f’^^r  time  histories  down  to  label  107,  and  for  TRIM  start  and  finish  tabula- 
tions dovni  to  label  109. 

The  variable  EH  is  filled  for  the  harmonic  einalysis  parameters  and  PDATA 
is  filled  for  the  time-history  parameters.  The  PDATA  being  collected  is 
either  for  the  normal  time  histories  being  plotted  during  the  trim  process 
or  for  the  time  histories  of  some  of  the  parameters  being  harmonically 
analyzed,  (uxtra  revolutions  are  added  to  trim  if  the  harmonic  analysis 
flag  is  on.)  The  definitions  of  some  key  variables  are:  TIMET,  the 

elapsed  time  in  trimj  PTII4E,  the  plot  time  same  as  TIMET  until  500  plot 
joints  are  collected;  NPTS,  the  *'umber  of  plot  points  collected;  NPLTf, 
the  nimiber  of  time  points  skipped  to  the  next  plot  pointy  and  IPTO,  a 
■cime  point  counter  between  plot  points.  SCALE  and  TSCLE  are  the  rotor 
revolutions  and  seconds  per  inch  of  plot  paper.  CYCLE  is  the  number  of 
revolutions  from  the  beginning  of  trim,  and  AZIMUTH  is  the  number  of  time 
points  f:om  zero  azimuth. 

After  some  fast  trim  code  from  label  109  to  108l,  the  program  conducts  some 
end  of  trii  tests.  The  first  test  is  simply  to  determine  if  CYCLE  equals 
TCuT,  the  total  number  of  cycles  the  progr  m is  allowed  to  trim,  an  input. 
All  ■the  other  tests  are  based  on  the  trim  variable  size.  The  trim  comple- 
tion tests  are  done  on-^'e  per  revolution,  and  permit  TRIM  to  exit  if  the 
change  in  each  and  every  control  (CONTRL)  from  one  time  point  to  another 
is  less  than  som."-  small  value  built  into  the  prograih.  The  trim  bomb  test 
causes  control  to  exit  TRIM  if  any  of  the  monitored  variables  exceed  a 
built-in  value.  This  portion  occupies  the  program  until  label  IlOx, which 
also  inc  odes  some  additional  fast  tri.m  code. 

Next  comes  thv:  pri.me  pa  t of  TRIM  where  the  trim  variables  are 
reevaluated.  Volume  III,  Sections  3.3.2  and  3.3.3  are  suggested  reading 
The  incremental  change  of  the  trim  variable  from  one  time  point  to  the  next 
jj  usually  based  on  accelerations  which  are  integrated  and  "filtered"  by 
the  FA  vai'iable.  In  effect,  the  high  frequencies  which  contain  the  vibra- 
tory components  are  smoothed  out  by  the  integration  process. 
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The  filter  also  has  bandpass  characteristics,  and  as  such  only  retains 
current  and  near  current  iteration  values.  Old  trials  are  washed  out.  In 
general  the  trim  variables  then  drive  the  integrated,  filtered  acceleration 
errors  to  zero.  For  the  flexible  shaft  option, the  shaft  angles  eire  deter- 
mined from  auxiliary  equations  giving  the  static  deflection  as  a function 
of  load.  The  static  deflections  thus  determined  are  integrated  and  fil- 
tered for  use  in  the  main  stream  trim  calf  lations. 

The  program  at  this  point  also  readjusts  the  pilot  stick  and  the  control 
gyro  positions  once  the  swashplate  trim  is  complete.  All  the  above 
occupies  the  program  down  to  label  117,  which  ends  the  TRIM  loop. 

Beyond  the  117  label  is  programming  whicn  reverses  the  usual  updating  of 
old  values  of  the  motions  of  the  primary  degrees  of  freedom.  The  purpose 
here  is  to  have  the  first  time  point  in  FLY  the  same  as  the  last  time  point 
in  TRIM.  After  label  lli*7  there  is  a short  portion  which  updates  and 
punches  new  trim  cards,  if  so  flagged,  and  calls  RCPLOT.  The  calls  for 
RCPLOT  between  labels  310  and  311  are  for  the  harmonic  analysis  plots.  The 
program  ends  by  computing  sine  and  cosine  components  as  well  as  the  magni- 
tude and  phase  of  certain  rotor  loads  and  the  feather  angle.  This  data  is 
tabulated  by  the  name  XYZ. 

1.3.31  TORSI 

The  isolated  twist  response  of  the  blade  at  every  station  is  determined  in 
this  subroutine.  TORSI  initializes  while  entry  TORS  does  the  twist  comput- 
ing at  every  time  point.  FX,  FY,  and  FZ  are  the  blade  shears  integrated 
from  the  tip  to  the  statiOii  in  question.  They  equal  the  root  values 
F(l4,I),  F(5,I),  and  F(6,l)  minus  the  values  from  FSV  saved  in  subroutine 
SWEEP  during  the  integration  proceeding  from  the  root  to  the  station  in 
question.  TORQUE,  then,  is  torque  at  the  shear  center  where  XSV,  YSV,  and 
ZSV  are  the  distance  to  the  shear  center  in  blade  root  axes  as  saved  in 
SWEEP.  DSOGJ  .is  the  reciprocal  of  the  torsion  stiffness,  DTHG  is  the 
increment  in  twist  angle  over  a blade  section, and  THTORS  is  the  total 
twist  from  the  tip  to  the  given  station.  Only  the  displacement  is  com- 
puted. The  t ist  velocity  THTRD  is  taken  as  zero  to  help  prevent  numerical 
instabilities  during  dynamic  stall  on  the  retreating  blade. 

3.3.32  TRMPUH 

TRMPUN  punches  those  inputs  which  may  vary  dixring  TRIM  and  describe  the 
trim  state  of  the  model  (Trim  Save  option).  These  variables  include  the 
trim  control  variabies  and  the  motions  of  the  degrees  of  freedom.  When 
these  trim  save  cares  are  added  to  the  input  deck,  the  next  case  starts  off 
where  the  old  one  ec.-^ed.  These  cards  can  shorten  the  time  to  trim  even  if 
new,  somewhat  diffeier.t,  flight  conditions  are  used. 
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& Entry  DSHIFT  updates  almost  the  same  set  of  inputs  as  TRMPUN.  The  update 

I is  saved  in  RAS.  As  such  DSHIFT  might  be  useful  during  a single  submittal 

^ of  a number  of  similar  runs  where  trim  save  card  are  not  available. 

I'  DSHIFT  and  TRMPUN  are  called  by  TRIM  if  RA(2000)  = TRMUPD  or  RA(i+T)  = 

I IPUNCH  flag  is  on  (=l). 

1.3.33  XTERMO 

i This  subprogram  contains  a major  portion  of  the  transformation  matrices 

required  in  SWEEP.  A given  entry  may  use  up  to  three  angles.  Multiple 
■ . Euler  rotations  as  well  as  derivatives  are  coded.  Individual  call  descrip- 

tions are  referred  to  the  description  in  SWEEP.  The  coding  of  this  sub- 
program has  been  machine  produced  by  the  IBiM  FORMAC  procedure. 

1.3.314  XTRPl,  XTRP2,  XTRP3 

This  is  a set  of  special  linear  interpolation  routines  designed  for  speed. 
The  interpolation  technique  does  not  require  the  classic  table  search  to 
find  the  data  of  interest.  These  routines  do  require  that  the  tabulated 
function  be  evaluated  at  constant  argument  intervals. 

Algoritim  - If  is  the  first  argument  valu*  and  the  data  is  evaluated  at 
Ax  intervals,  then  the  location  within  the  table  of  bracketing  function 
values  for  a value  of  the  argument  x is  n where 


n = 


(1~93) 


The  brackets  indicate  the  largest  integer  whose  magnitude  does  not  exceed 
the  magnitude  of  (x  - Xq)/  Ax.  The  interpolating  formula,  based  on  the 
point-slope  formula,  is  given  as 


y 


y + 

n 


dx(y 


n+1 


(I-9M 


where 


(1-95) 
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This  is  shown  to  he  correct  if 


x-x.  = ri*Ax  + (x-x) 
0 n 


(1-96) 


is  substituted  in  the  equation  for  dx. 

The  quantities  Xq  and  Ax  are  stored  with  the  table  of  function  values. 
The  function  of  each  routine  is 


XTRPl 

computes  F = f(x) 

XTRP2 

computes  F = f(x,y) 

XTRP3 

computes  F = f(x,y,z) 

There  is  one 
contain  the 

general  table  format  for  all  routines.  The  table,  T,  must 
following : 

T(l) 

^0 

T(2) 

number  of  x breakpoints  - 1 

T(3) 

1/Ax 

T(h) 

^0 

T(5) 

number  of  y breakpoints  - 1 

T(6) 

1/^y 

T(7) 

^0 

T(8) 

number  of  z breakpoints  - 1 

T(9) 

1/Az 

T(10),. 

..  table  data  on  increasing  x,  then  y,  and  then  z. 

1.3.35  XTRP 

XTRP  is  a generalized  table  driven  interpolation  routine  capable  of  evalu- 
ating functions  of  1,  2,  or  3 variables,  either  linearly  or  quadratically. 
The  use  of  XTRP  is  thoroughly  explained  with  comment  cards  associated  with 
the  module. 

1.3.36  XTRPl4  and  XTRPCA 

Together,  these  routines  compute  main  rotor  blade  aerodynamic  coefficients 
C,  and  Cjj  where 
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Cl  = Cl(M,  ot.  Cl_,  t/c) 


% = %(M,  <x,  , t/c) 

i 


These  functions  are  defined 


over  the  ranges: 


C S a < 36o 
0.1  s M < .9 


0 i Cl  < .69 

i 


0.06  S t/c  < , ■ 


interpola- 

into  the  software  via  block  data  subprograms 
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2.  COMON/SUBROUTINE  DIRECTORY 


The  collection  of  subprograms  which  constitute  REXOR  contain  many  cross- 
referenced  COMMON  blocks  and  subroutines.  Table  2-1  is  presented  as  an 
aid  to  developing  the  source  and  usage  of  any  particular  item.  The  verti- 
cal listing  on  this  table  gives  all  the  routines,  and  names  or  unnamed 
(one)  COMMON  blocks  used  in  REXOR.  A second  column,  headed  by  'T'  gives 
the  use  of  the  entry.  The  coding  used  is  *M'  for  main  program,  'S'  for 
subroutine,  'E'  for  entry  of  a subroutine,  'C  for  common  block.  The  list 
of  names  is  in  alphabetical  order  and  numbered.  The  alphabetizing  is  by 
main  to  subroutines  with  a subalphabetizing  of  entry  points  under  each 
subroutine.  The  COMMON  blocks  are  listed  last  except  for  subroutines  not 
included  in  the  source  deck.  These  subroutines  are  usually  part  of  the 
computer  operating  (precompiled)  package,  and  not  particularly  associated 
with  REXOR.  The  numbering  is  repeated  horizontally,  and  corresponds  to 
the  vertical  name  list. 

The  vertical  list  on  the  left-hand  side  is  the  calling  or  active  routine  or 
element,  and  the  horizontal  line  lists  the  routines  called  or  referenced. 
Numbers  at  grid  intersections  show  there  is  a reference  and  the  level  of 
reference.  One  indicates  a direct  reference.  Two  or  three  show  there  is 
one  or  two  intermediate  references,  respectively.  Note  that  a subroutine 
name  will  show  all  the  references  to  all  the  entry  points  bounded  by  that 
name.  By  elimination,  the  references  associated  with  the  subroutine  romc 
up  to  the  point  of  the  first  entry  nasne  can  be  determined. 
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TABLE  _-l.  COMMON  SUBROUTINE  DIRECTORY  - Continued 


3.  TIME  MD  SPATIAL  INTEGRATION 


The  equations  of  Volume  I have  treated  the  flying  vehicle  to  be  composed 
mostly  of  a number  of  rigid  body  elements.  For  the  blade,  however,  con- 
venience dictates  an  integration  process  with  distributed  mass  and  inertial 
data.  The  subroutine  SWEEP  ends  with  a trapezoidal  integration  of  all 
blade  generalized  masses  and  loads.  This  algorithm  approximates  th  inte- 
gral from  root  to  a station  X. 


F(X) 


X 


/ 


F'  dX 


(2-1) 


by 


i=l 


[F*(i)  + F»(i  + 1)]  [X(i  -*•  1)  - X(i)] 
2 


(2-2) 


The  more  complex  Simpson  ruje  is  frequently  used,  but  trapezoida]  integra- 
tion is  more  convenient  because  of  its  ready  geometric  interpretation,  and 
is  especially  useful  for  applying  aerodynamic  'tip  loss'. 

The  time  integration  of  the  degrees  of  freedom  stai’ts  with  a Euler  inte- 
gration which  simulates 


as 


P(t) 


P dt  + P(0) 


(2-3) 


P(nt) 


p(j)  At  + p(o) 

j=i 


(2-M 
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At  being  a constant  interval  between  time  points.  After  four  time 
points  are  computed,  and  enough  old  data  collected,  a switch  is  made  to 
Adams-Bashforth  integration. 


P(Aj) 


"t 

[55P(j)  - 59P(J  - 1)  + 37P(J  - 2)  - 9P(J  " 3)  ] + P(0) 

J=1 

(2-5) 


This  sequence  is  used  in  INTO  for  the  majority  of  the  degrees  of  freedom, 
and  in  BIRDl  for  gyro  and  shaft  bending  motions. 


Note  that  the  operation  of  the  integration  scheme  can  severely  affect  the 
apparent  mode  damping.  This  is  important  to  the  inplane  mode  which  char- 
acteristically has  low  damping.  The  question  as  to  whether  the  mode  is 
stable  or  unstable  is  of  extreme  interest.  For  example,  suppose  computa- 
tions are  180  points  per  revolution,  or  two  degree  azimuth  in  FLY.  This 
suggests  that  the  phase  of  a parameter  oscillating  at  one  per  revolution 
(IP)  is  limited  to  an  ultimate  resolution  of  two  degrees;  at  2P,  ultimate 
resolution  is  four  degrees,  etc.  For  a second -order  system  undergoing 
transient  response 


X + 


gw^x 


+ CD, 


= f 


(2-6) 


the  ratio  of  the  magnitude  of  the  damping  term  to  either  the  spring  or 
mass  term  for  an  oscillation  at  the  natural  frequency  uq  is  g.  In  the 
example  at  IP,  g would  be  limited  in  precision  to  about  0.01  (-  2 deg). 
This  level  is  roughly  the  order  of  structural  damping.  The  analysis  of  a 
large  number  of  oscillations  might  give  greater  precision;  however,  the 
above  position  is  conservative.  The  foregoing  argument  assumes  the  model 
has  all  the  physical  components  of  importance  described  and  a sufficient 
number  of  blade  mass  elements. 
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k.  AERO  LOOKUP  - MAIN  ROTOR 

The  sequence  of  main  rotor  blade  section  aerodynamics  computations  is 
explained  in  Volume  I Sections  7.2.2  through  T.2.U.  The  computer  code 
corresponding  to  these  procedures  is  explained  in  this  section. 

All  of  these  computations  are  initiated  through  the  routine  SWEEP. 

Depending  on  the  selection  of  ILOOK  = RA(2689),  either  the  direct  aero 
tables  (seven  table  lookup)  or  tailored  to  specific  geometry  tables  (fast 
aero)  are  used.  Note  that  these  tables  can  be  replaced  by  user  supplied 
aerodynamic  data  so  long  as  the  CALL  arguments  coincide. 

. ].  BASE  SET  OF  AIRFOILS.  CAMBER  AND  THICKNESS  RATIO 

The  seven  table  lookup.  Table  1*-1  is  a grouping  of  the  available  table 
set,  Table  h-2,  according  to  thickness  ratio  and  design  lift  coefficient. 
SWEEP  accesses  the  Cj^,  Cp  part  of  this  data  for  ILOOK  = 1 via  the  inter- 
polation routine  XTRPlj.  XTRPI4  logic  determines  the  proper  design  lift 
coefficient  and  thickness  tables  to  bracket  the  desired  values  through  the 
routine  XTRPCA,  and  performs  the  interpolation.  The  data  is  sorted  and 
retrieved  by  the  COMMON  naming.  Table  J»-2. 

The  Cj^i  value  required  is  determined  through  the  routine  CMLOOK.  No  inter- 
polation to  design  Cp  or  thickness  ratio  is  made.  Rather,  the  input  IFOIL  = 
RA(2690)  selects  0012  airfoil  (=0)  or  23012  (=l). 
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k.2  REDUCTION  TO  SMALLER  TABLE  SETS 

The  seven  table  data  (plus  C^  tables)  may  be  collapsed  to  a set  directly 
keyed  to  the  geometry  at  hand.  The  scheme  used  in  REXOR  uses  equally 
spaced  data  index  points,  which  leads  to  rapid  access.  This  procedure, 
called  Fast  Aero,  is  called  from  SWEEP  for  ILOOK  = 0.  SWEEP  uses  the  Fast 


Aero  tables  through  the  routine  AERO.  This  routine  is  also  used  by  STALL 
(dynamic  stall)  which  is  called  by  SWEEP  for  ISTALL  = RA(2555)  = 0. 

The  Cp  and  Cp  data  for  Fast  Aero  is  grouped  according  to  normal  angle  of 
attack  range  (j^30  degress)  with  Mach  effects  and  the  remaining  angle  dis- 
tribution treated  as  incompressible  flow.  The  Fast  Aero  example  is  shown 
in  Figure  ^*-1,  and  is  constructed  and  referenced  by  thickness  ratio  for 


TABLE  U-1.  SEVEN  TABLE  LOOKUP 


(: 


Table 

Grouping 

Included 

Tables 

Design 

CL 

Thickness 

1 

CMnCL 

CMTICD 

CT3SCL 

CT3SCD 

0 

0.12 

2 

CMT2CL 

CMT2CD 

CT2SCL 

CT2SCD 

0 

0.09 

3 

CMT3CL 

CMT3CD 

CTISCL 

CTISCD 

0 

0.06 

1* 

CMTltCL 

CMTitCD 

CTl*SCL 

CTl*SCD 

0.09 

0.06 

5 

CMT5CL 

CM1’5CD 

CT5SCL 

CT5SOD 

0.39 

0.09 

6 

cmt6cl 

cmt6c.o 

ct6scl 

ct6scd 

0.35 

0.12 

7 

CMT7CL 

CWr7CD 

CT7SCL 

CTiSCD 

0.69 

0.12 

i 
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TABLE  U-2.  BASE  AIRFOIL  DATA 


Table 

Common 

Name 

Camber 

Thickness 

Design 

CL 

a 

Range 

Mach 

Range 

Output 

NACA  Type 

1 

CMTICL 

0 

0.12 

0 

+30 

(1) 

CL 

(0)3012 

CMTICD 

0 

0.12 

0 

±30 

(U) 

CD 

(0)3012 

2 

CMT2CL 

0 

0.09 

0 

±30 

(2) 

CL 

(0)3009 

CMT2CD 

0 

0.09 

0 

+30 

(3) 

CD 

(0)3009 

3 

CMT3CL 

0 

0.06 

0 

+^30 

(5) 

CL 

(0)3006 

CMT3CD 

0 

0.06 

0 

+30 

(6) 

CD 

(0)3006 

1; 

CMTUCL 

0.6 

0.06 

0.09 

+30 

(5) 

CL 

(0.6)3006 

CMTUCD 

0.6 

0.06 

0.09 

+30 

(6) 

CD 

(0.6)3006 

5 

CMT5CL 

2.6 

0.09 

0.39 

±30 

(T) 

CL 

(2.6)3009 

CMT5CD 

2.6 

0.09 

0.39 

+ 30 

(8) 

CD 

(2.6)3009 

6 

cmt6cl 

2.3 

0.12 

0.35 

+30 

(D 

CL 

(2.3)3012 

cmt6cd 

2.3 

0.12 

0.35 

±30 

(9) 

CD 

(2.3)3012 

T 

CMTTCL 

h.6 

0.12 

0.69 

+30 

(7) 

CL 

(i*. 6)3012 

CMT7CD 

k.e 

0.12 

0.69 

±30 

(9) 

CD 

(1*. 6)3012 

8 

CTISCL 

0 

0.06 

0 

+30 

+330 

INCOMP 

CL 

(0)3006 

CTISCD 

0 

0.06 

0 

+30 

+330 

INCOMP 

CD 

(0)3006 

CT2SCI. 

0 

0.09 

0 

+30 

+330 

INCOMP 

CL 

(0)3009 

CT2SCD 

0 

0.09 

0 

+30 

+330 

INCOMP 

CD 

(0)3009 

CT3SCL 

0 

0.12 

0 

+30 

+330 

INCOMP 

CL 

(0)3012 

CT3SCD 

0 

0.12 

0 

+30 

+330 

INCOMP 

CD 

(0)3012 

9 

CTi*SCL 

0.6 

0.06 

0.09 

+30 

+330 

INCOMP 

CL 

(0.6)3006 

CTitSCD 

0.6 

0.06 

0.09 

+30 

+330 

INCOMP 

CD 

(0.6)3006 
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TABLE  h-2  - Continued 


Table 

Common  Design 

Name  Camber  Thickness  CL 

0 

Range 

Mach 

Range  Output 

NACA  Type 

9 

CT5SCL  2.6 

0.06 

0.39 

+30 

+330 

INCOMP 

CL 

(2.6)3006 

CT5SCD  2.6 

0.06 

0.39 

+30 

+330 

INCOMP 

CD 

(2.6)3006 

CT6SCL  2.3 

0.12 

0.35 

+30 

+330 

INCOMP 

CL 

(2.3)3012 

CT6SCD  2.3 

0.12 

0.35 

+30 

+330 

INCOMP 

CD 

(2.3)3012 

CTTSCL  ii . 6 

0.12 

0.69 

+30 

+330 

INCOMP 

CL 

(1.6)3012 

CTTSCD  U.6 

0.12 

0.69 

+30 

+330 

INCOMP 

CD 

(1.6)3012 

10 

CMTICM  2 

0.08 

0.3 

+30 

(10) 

CM 

23008 

CTISCM*  0 

0.12 

0.0 

+30 

+330 

INCOMP 

CM 

(0)3012 

11 

CMT2CM  0 

0.12 

0.0 

±30 

(11) 

CM 

(0)3012 

CTISCM*  0 

0.12 

0.0 

+30 

+330 

INCOMP 

CM 

(0)3012 

^REPEATED 

Mach  Ranges 

(1)  0,  0.1,  0.2, 
0.85,  1.0 

0.3, 

0.1,  0.5, 

0.55, 

0.6,  0.65, 

0.7,  0. 

.75, 

(2)  0,  0.1,  0.2, 
0.35,  1.0 

0.3, 

C.l,  0.5, 

0.55, 

0.6,  0.65, 

0.7,  0.757, 

(3)  0,  0.3,  O.ii, 

0.5, 

0.6,  0.7, 

0.8,  0.915,  1 

(k)  0,  0.738,  0.78,  0.805,  0.85 

, 0.9 

(5)  0,  0.2,  O.li, 

0.5, 

0.6,  0.65 

, 0.7, 

0.75,  0.805 

i,  0.85, 

, 1 

(6)  0,  0.1,  0.5, 

0.6, 

0.7,  0.75 

, 0.8, 

0.85,  0.9, 

0.935, 
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TA.''LE  h~2  -■  Concluded 

Mach  Ranges 

(7) 

0,  0.2,  0.3, 

O.h, 

0.5,  0.55,  0.6,  0.6f 

, 0,7,  0.757,  0.85,  i.O 

(8) 

0,  0.3,  0.l4, 

0.5, 

0.6,  0.7,  0.8,  0.85, 

0.915,  1 

(9) 

0,  0.2,  0.I4, 

0.5, 

0.6,  0 .7,  0.8,  0.9 

(10) 

0,  0.3,  0.I4, 

0.925,  1 

0.5, 

0.55,  0.6,  0.65,  0.7 

, 0.75,  0.8.  0.875, 

(11) 

0,  0.1,  0.2, 
0.8,  0.9 

0.3, 

0.I4,  0.5,  0.55,  0.6, 

0.65,  0.7,  0.75, 

Due  to  the  limited  section  data  available,  the  data  is  not  displayed  as 
a function  of  thickness  ratio  in  Fast  Aero,  but  is  bivariantly  interpolated 
for  +30  degrees  and  univariantly  interpolated  for  incompressible  flov, 
XTRPl.  The  base  data  is  for  a design  lift  coefficient  of  0.3  and  thick- 
ness ratio  of  0.08  (NACA  23008).  The  compressible  flow  range  part  of  this 
data  is  curve  shifted  on  the  angle  of  attack  axis  to  account  for  other 
thickness  ratios  as  shown  in  Figure  It-C.  This  operation  is  performed  in 
AERO. 


k.3  RANGE  OF  ENTRIES 

The  existing  Fast  Aero  set  is  specifically  constructed  for  the  AH-56A 
blade.  However,  the  seven  table  set  contains  a large  spectrum  of  sections 
made  Ijy  drooping  the  nose  portion  of  symmetrical  four  digit  foils.  These 
are  known  as  the  five  digit  series  and  the  equivalent  designation  is  given 
on  Table  h-2.  Therefore  a Fast  Aero  set  could  be  constructed  for  any 
blade  composed  of  five-digit  airfoil  sections. 


k.k  DATA  CORRECTIONS 

A Cm  trim  is  provided  in  AERO  in  the  comp  essible  flow  range  for 
(IBLADE,  liA  1300)  = 3 by  the  coefficient  DCMRl. 

/.  second  trim  is  available  in  SWEEP  to  account  for  blade  trim  tab. 
This  trim  quantity  DCMR  operates  between  blade  stations  KTO  and  KTl . 

An  override  of  Cj^  = 0,  Cm  = 0 is  set  in  SWEEP  for  the  blade  root  cutout. 
The  Cp  value  calculated  is  used  as  calculated  from  the  tables. 
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COMMON  ANGLE  OF  ATTACK  THICKNESS  RATIO  MACH 


TABLE 

NAME 

OlTTPLnr  MIN 

MAX 

A 

MIN 

MAX 

A 

MIN 

M-Ar 

A 

1 

CC-ATM 

1 

CL 

-36.6 

33.4 

1.0 

0.08 

0 12 

0.03 

0 

0.9 

0.05 

2 

CLAT 

CL 

33.33 

3m.37 

0.5 

0.06 

0.12 

0.03 

0 

0 

1 

3 

CDATM 

CO 

-me 

33.4 

1.0 

0.06 

0.12 

0.03 

0 

0.9 

0.05 

4 

COAT 

CD 

33.38 

330.37 

0.5 

0.06 

0.12 

0.03 

0 

0 

1 

5 

CMAM 

C»l 

-30 

30 

0.5 

N/A 

N/A 

N/A 

1 

0.3 

1 

0.9 

0.025 

8 

CMA 

CM 

30 

330 

0.5 

N/A 

N/A 

N/A  1 

1 0 

0 

1 

NOTE  1;  FOR  CL.  CO  DATA  ACCESSION  IS  BY  MEANS  OF  THICKNESS  RATIO  AS  THR.  DESIGN 
CL  AND  BLADE  STATIOfi  ARE  LINEARLY  RELATED 


CAMBER 

o.e^ 


DE3 


4.6 


NOTE  2:  CM  DATA  IS  TABULATED  TO  ANGLE  OF  ATTACK  AND  MACH  (OPT)  ONLY.  DATA  IS 
CURVE  SHI  STED  FOR  THICKNESS  RATIO  {!.£.,  STATION)  FUNCTION 


Figure  !*-!.  Fast  Aero  Tables. 
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SECTION  MOMENT  COEFFICIENT 


SECTION  MOMENT  COEFFICIENT  VS  ANGLE 
OF  ATTACK  {-30®  TO  30«)  FOR  ALL  SECTION 
LOCKHEED  CHORD  REFERENCES 


0.926  TO  1.0  OTHER  SECTIONS.  PHASE  til  SLADES 

0.875  ‘^EAD  ■ ® ^ *'^2  *^3  ® 

0.80  WHERE 

0.76  K^  - 243,  t/c  > 0.03 

S.  - 400,  t/c  < 0.08 

Kj-I.  |fll<10° 

. ^ - 130- |a|l/14,  16O<|oK30O 

N.  \ Kj  * <M  - 0.7)  /0.7,  M < 0.7 

V X.  \ - 0,  M > 0.7 

. X \ \ ACCOUNTS  FOR  t/c  AND  C,  DISTRIBUTION 
^ X \\  ‘-DES 


C|^  ® C|^  "t"  0.01 

^ ""READ 


0 TO  0.30 


0 TO  0,30 


NACA  23008 
BASIC  DATA 

t/c  “ 0.08 


r/R  - 0.67 


,0.925  TO  1.0 


ANGLE  OF  ATTACK.  Qread  " 


Figure  h-2.  Fast  Aero  Calculations. 
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5.  COMPLgTE  SOURCE  LISTING 


Due  to  the  iai'ge  number  of  pages  in  the  REXOR  program  source  listing,  this 
material  is  handled  under  a separate  binder.  Copies  may  be  obtained  from 
the  distributing  agency,  USAAMRDL,  Eustis  Directorate,  Ft.  Eustis, 
Virginia. 
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LIST  OF  SYMBOLS 


SYMBOL, S 


a 

^0 

•‘^1 

[a] 

^.2.3 

^2n 


IS 


IS 


[c3 


'M 


arbitrary  vector 

2 

acceleration  vector,  ft/sec 

longitudinal  component  of  blade  first  harmonic  flapping,  rad 
generalized  mass  element  matrix 
modal  variables 

generalized  displacement  of  n^  blade,  first  mode 
generalized  displacement  of  nth  blade,  second  mode 
generalized  displacement  of  n^  blade,  third  mode 
cosine  component  of  blade  first  harmonic  cyclic,  rad 
number  of  main  rotor  blades;  arbitrary  vector 
dissipation  function 

sine  component  of  blade  first  harmonic  cyclic,  rad 
blade  segment  chord,  ft 
damping  matrix 

aerodynamic  drag  coefficient 
aerodynamic  lift  coefficient 
aerodynamic  pitching  moment  coefficient 
power  coefficient 
thrust  coefficient 
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^X,Y,Z 

C(k) 

d 

dr 

dt 

a/dt 

(d/e)o 

(d/e)^ 

e 

El 

^'iMR 

F 


f' 

x.y.z 

F 

(jj  ,0 

F 

6PH 

g 

®X,Y,Z 

G 

{G} 

G 

GJ 


lirjfjar  damping,  lb/t‘t/sec 

rotar7  damping,  ft-lb/rad/sec 

blads  bending  to  feathering  couplings 

lift  deficiency  function 

infinitesimal  increment 

increment  in  rotor,  radius,  ft 

increment  in  time,  sec 

derivative  with  respect  to  time 

swashr’. ate  to  feather  gear  ratio,  zero  collective 

swashplate  to  feather  gear  ratio  slope  with  collective 

pitch  horn  effective  crark  arm,  ft 

blade  bending  stiffness  distribution,  Ib-ft^ 

ground  effect  factor  for  main  rotor 

factor;  force,  ib 

force  components  along  X,Y,Z  directions,  lb 

generalized  force  about  (p , 0 , axis 

feathering  mode  generalized  force 
2 

gravity,  ft/sec 

gravity  components  along  X,Y,Z  directions 
gear  ratio 

generalized  force  vector 

gyro  angular  acceleration  partial  product 

2 

blade  torsional  stiffness,  Ib-ft 

= Em.  X.^,  I’lug-ft^ 

11 

= Em.  Y.^,  slug-ft^ 

11 
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'z 

= 

Z^^,  slug-ft^ 

^xx 

= 

Zi  ),  slug- ft' 

Iyy 

= 

Zm. 

1 

Zi^),  slug-ft 

^zz 

= 

Zm^ 

(Xi^. 

Yi^"),  slug-ft 

T 

“Xi' 

= 

Zm. 

1 

Xi  Yi, 

slug-ft^ 

^xz 

= 

z^ 

Xi  Zi, 

2 

slug-ft 

'yz 

= 

Zm. 

1 

Y.  Z., 

slug-ft^ 

i unit  vector 

J unit  vector 

J advance  ratio 


X.Y.Z 


numJ’jer  of  blade  radial  stations;  reduced  frequency, 
rad/sec;  unit  vector 

spring  matrix 

blade  spring  matrix  element 

spring  constants  along  X,Y,Z  direction,  Ib/ft 

spring  rates  about  d,  li'  axis,  ft-lb/rad 

location  inboard  feather  bearing,  ft 

location  outboard  feather  bearing,  ft 

radial  location  of  intersection  of  precone  and  feather 
axis,  ft 

tension  torsion  pack  length,  ft 

rolling  moment,  ft-lb 

mass  of  element,  slugs 

summed  fuselage  coordinate  mass , slugs 

summed  hub  axis  ii3ss , slugs 

mass  of  ith  particle  or  blade  segment,  slugs 
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NA 


t 


T 

[t] 


u 

”c,P,S,T 

V 


w 


''iMR 

u 

iTR 

X 

X 


sw 


TTO 


1,2,3 


ONA 


Laplace  variable,  path  of  motion  of  particle  p 
blade  spline  length  along  neutral  axis  locii,  ft 
time 

kinetic  energy,  ft-lb 
transformation  of  coordinates  matrix 
tension  in  tension  - torsion  pack,  lb 
velocity  in  X direction,  ft/sec 

pi:tential  energj'  function,  ft-ib;  strain  energy,  ft-lb 

air  velocity  on  blade  element,  ft/fc 

velocity  in  Y direction,  ft/sec 

trajectory  velocity 

velocity  in  Z direction,  ft/sec 

main  rotor  collective  inflow,  ft /sec 

tail  rotor  collective  inflow,  ft/sec 

motion  in  X direction,  ft;  blade  span  location 

coordinate  direction;  axis;  deflection,  ft;  location,  ft; 
cross  product 

blade  radial  station  oi'  sweep  and  Jog,  ft 
trajectory  path,  ft 
tail  rotor  longitudinal  force,  lb 
motion  in  Y direction,  ft 

coordinate  direction;  axis;  deflection,  ft;  location,  ft 
tension  torsion  pack  outboard  end  modal  coefficients 

difference  between  Y direction  locations  of  eg  and  neutral 
axis  points  of  blade  element,  ft 
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■“sp 

M 

[M] 

«rk 

^.Y,Z 

N 


P 

^iMR 

q 


^iMR 

Q 

QLOADS 

S’R 


R 

R. 


Z(Ji,Ze 


swashplate  summed  mass , slugs 

pitching  moment,  ft-lb;  = Em^,  slugs 

generalized  mass  matrix 

generalized  mass  matrix  element 

= Era.  X. , slug-ft 
11 

= Em.  Y.,  slug-ft 
11 

= Era.  Z.,  slug-ft 
11 

moments  about  X,Y,Z  axis,  ft-lb 
blade  torsional  moment,  ft-lb/ft 
number  of  system  particles 

angular  velocity  about  X axis,  rad/sec;  particle 

main  rotor  pitch  moment  inflow,  ft/sec 

generalized  coordinate;  angular  velocity  about  Y axis, 
rad/sec 

main  rotor  roll  moment  inflow,  ft/sec 
generalized  forcing  function 

aerodynamic  pressure  times  reference  wing  area,  lb 
total  nonmain  rotor  aerodynamic  loads  matrix 
tail  rotor  torque,  ft-lb 

general  vector;  radius  of  curvature,  ft;  angular  velocity 
about  Z axis,  rad/sec;  notation  for  (X,Y,Z) 

static  blade  shape 

vector  displacement  of  particle  p in  X,Y,Z  axis  system 
vector  displacement  of  x,y,z  origin  in  X,Y,Z  system 
gyro  damper  coupling  ratios 
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SP 


TTO 


1,2,3 


^OBL 


OF 

^OSP 


"OTTI 


FA 


PHn 


3TR 

3c/9o 


A 

P 


motion  in  Z direction 

coordinate  direction;  axis;  deflection,  ft;  location,  ft 

relative  swashplate  vertical  displacement  with  respect  to 
the  hub,  ft 

tension-torsion  pack  outboard  end  modal  coefficients 

teetering  rotor  undersling,  ft 

hub  set  distance  above  fuselage  set,  ft 

hub  set  distance  above  swashplate  set,  ft 

blade  vertical  offset  at  outboard  end  of  tension  - torsion 
pack,  ft 

angle  of  attack,  rad 

angle  of  attack  with  hub  set,  rad 

sideslip  angle,  rad 

blade  feathering  angle,  rad 

feathering/pitch-horn  bending  or  dynamic  tc.rsion 
generalized  coordinate  displacement 

blade  droop  relative  to  precone  angle,  rad 

blade  sweep  angle,  rad;  dynamic  stall  delay,  sec 

trajectory  path  angle  with  E set,  rad 

limit  deflection,  rad;  freeplay,  rad;  small  increment 

tail  rotor  pitch  - flap  coupling 

dovnwash  factor  of  wing  on  horizontal  tail 

vector  notation  of  4>,  0,  \l/ 

rotation  about  Y axis,  rad 

collective  blade  angle,  rad 

sideslip  at  blade  element,  rad 

3 

ai.r  density,  slugs/ft 


88 


T 


T 


0 


"W 


3 

SUBSCRIPTS 


time  constant,  sec;  natureil  period,  sec 
feathering  axis  precone,  rad 
rotation  about  X axis,  rad 
feathering  angle,  rad 

feathering  angle  of  blade  element  of  nth  blade,  rad 

blade  root  reference  feather  angle,  rad 

blade  torsion,  rad 

sum  of  blade  twist  and  torsion,  rad 

wake  angle  of  main  rotor,  deg 

rotation  about  Z axis,  rad;  sideslip  angle  with  hub  set,  rad 

control  input  axis  rotation  from  swashplate,  rad 

pitch  lead  angle,  deg 

trajectory  path  yaw  with  E set,  rad 

main  rotor  apparent  airflow  angle,  rad 

rotational  speed,  rad/sec;  angular  velocity,  rad/sec; 
natural  frequency,  rad/sec 

partial  derivative,  derivation 


a arbitrary  coordinate  set  a 

A due  to  aerodynamics 

b arbitrary  coordinate  set  b 

BEND  associated  with  blade  elastic  bending 

BLE  blade  element  coordinate  system 

BLn  blade  reference  axis  system  for  the  nth  blade 
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c 

CG 

CORK 

DW 

DYN 

E 

ENG 

EST 

F 

FA 

FB 

Fn 

FR 

G 

GEN 

GFB 

GSP 

GUB 

H 

HT 

i 

IB 

J 

J 


associated  with  pilot  control  input,  chordwise 

associated  with  center  of  gravity  location 

corrective,  correction 

referring  to  downwash 

referring  to  dynar.ic  component 

earth  axis 

associated  with  powerplant  - engine 
estimated 

fuselage  axis;  associated  with  blade  feathering 

referring  to  blade  feather  axis 

associated  with  feedback 

associated  with  feathering  of  the  nth  blade 

due  to  friction 

referring  to  gyro  or  gyro  coordinate  system 

associated  with  gas  generator  section  of  powerplant 

associated  with  gyro  control  feedback 

gyro  to  swashplate  connection 

relating  to  gyro  gimbal  unbalance 

referring  to  hub  or  principal  reference  axis  system 

associated  with  horizontal  tail 

referring  tc  inflow,  particle 

referring  to  inboard  feather  bearing  location 

spring  matrix  index 

associated  with  blade  attachment  Joggle 
associated  with  gyro  end  of  feedback  rod  linkage 
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Jn 

k 

LAG 

LIMIT 

m 

MR 

n 

NA 

NEW 

NO 

NR 

OB 

OLD 

p 

PH 

r 

R 

REF 

RM 

S 

SC 

SP 

SP 

c 

S,  SP 


associated  with  feedback  rod  coming  from  the  nth  blade 

generalized  mass  index 

associated  with  lead-.’ag  damper 

signifying  limiting  value 

blade  mode  index,  spring  matrix  index 

associated  with  main  rotor 

blade  number  index 

referring  to  blade  segment  neutral  axis 

newly  determined  value 

normal  (to  airflow)  component 

pertaining  to  nonrotating  value 

referring  to  outboard  feather  bearing  location 

value  from  previous  time  step 

associated  with  propeller;  perpendicular  blade  component 

referring  to  pitch  horn 

generalized  mass  index 

referring  to  rotor  axis  system 

associated  wit'  ''lade  feather  reference  value 

referring  to  control  gyro  feedback  lever  moment 

referring  to  blade  spanwise  velocity;  general  mode;  static; 
structural;  shaft 

referring  to  blade  segment  shear  center 

referring  to  swashplate 

command  to  swashplate 

referring  to  swashplate  limit  stop 
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STEADY  steady  component 

SW  referring  to  blade  sweep  angle  location 

T associated  with  trajectory  path  relating  to  E axis; 

tangential  blade  component;  blade  torsion;  blade  twist 

TR  associated  with  the  tail  rotor 

TRIM  initial  or  trim  value 

TT  associated  with  tension  torsion  pack 

TTl  referring  to  inboard  end  of  tension  torsion  pack 

TTO  referring  to  outboard  end  of  tension  torsion  pack 

TW  associated  with  blade  twist  (built  in) 

UB  relating  to  control  gyro  unbalance 

UNSTEADY  associated  with  unsteady  component 

VT  associated  with  vertical  tail 

WING  associated  with  the  wing 

X relating  to  component  in  X direction 

Y relating  to  component  in  Y direction 

YA  relating  to  aerodynamic  component  in  Y direction 

Z relating  to  component  in  Z direction 

ZA  relating  to  aerodynamic  component  in  Y direction 

0 f nought)  associated  with  collective  value,  coordinate  axis 

value,  with  respect  to  principed  reference  eixis,  blade 
root  summation 


1,^,3 

with  respect  to 

IS 

first  harmonic  ' 

c 

with  respect  to 

3/i*  c 

with  respect  to 

92 


associated  the  feathering  mode  of  the  ntt  blade 
relating  to  component  in  the  ^ direction 
relating  to  component  in  the  6 direction 
relating  to  component  in  the  direction 

referring  to  inertial  reference 

matrix  transpose 

(bar)  average  quantity 

(prime)  slope  with  respect  to  blade  span 

(dot)  time  derivative  of  basic  quantity 

(double  dot)  second  time  derivative 

natrix  inverse 

vector  quantity 


blade  radial  station  index 
blade  number  index 
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